Results Replication
Neel Roy
2023-06-22
Here is my code:
Reformatting/Cleaning Up Data
Raw Counts File
cnts <- read.csv('raw_counts_Cooper_data_final.csv')
rownames(cnts) <- cnts$X
cnts$X <- NULL
cnts <- as.matrix(cnts)This takes raw_counts file and corrects the columns and index of the dataframe, then converts the file from a data frame to a matrix.
Meta Data
column_data <- read.csv('Cooper_Data_Meta.csv')
rownames(column_data) <- column_data$Sample
column_data <- column_data[,c("Tissue","Treatment","Stage")]
#column_data$Tissue <- factor(column_data$Tissue)
column_data$Treatment <- factor(column_data$Treatment)
#column_data$Stage <- factor(column_data$Stage)
#column_data$Treatment <- NULL
#column_data$Stage <- NULL
#column_data$Sample <- NULL
column_data## Tissue Treatment Stage
## SRR10912052 Serous_EOC Control 4B
## SRR10912053 Serous_EOC Platinum 4B
## SRR10912054 Serous_EOC Control 4
## SRR10912055 Serous_EOC Platinum 4
## SRR10912056 Serous_EOC Control 3C
## SRR10912057 Serous_EOC Platinum 3C
## SRR10912058 Serous_EOC Control 3C
## SRR10912059 Serous_EOC Platinum 3C
## SRR10912060 Serous_EOC Control 3C
## SRR10912061 Serous_EOC Platinum 3C
## SRR10912062 Serous_EOC Control 4
## SRR10912063 Serous_EOC Platinum 4
## SRR10912064 Serous_EOC Control 3C
## SRR10912065 Serous_EOC Platinum 3C
## SRR10912066 Serous_EOC Control 3C
## SRR10912067 Serous_EOC Platinum 3C
## SRR10912068 Serous_EOC Control 3C
## SRR10912069 Serous_EOC Platinum 3C
## SRR10912070 Serous_EOC Control 4
## SRR10912071 Serous_EOC Platinum 4
## SRR10912072 Serous_EOC Control 3C
## SRR10912073 Serous_EOC Control 3C
## SRR10912074 Serous_EOC Platinum 3C
## SRR10912075 Serous_EOC Control 3C
## SRR10912076 Serous_EOC Platinum 3C
## SRR10912077 Serous_EOC Platinum 3C
## SRR10912078 Serous_EOC Control 3C
## SRR10912079 Serous_EOC Platinum 3C
## SRR10912080 Serous_EOC Control 3C
## SRR10912081 Serous_EOC Platinum 3C
## SRR10912082 Serous_EOC Control 4
## SRR10912083 Serous_EOC Platinum 4
## SRR10912084 Serous_EOC Control 3C
## SRR10912085 Serous_EOC Platinum 3C
## SRR10912086 Serous_EOC Control 3C
## SRR10912087 Serous_EOC Platinum 3C
## SRR10912088 Serous_EOC Control 3C
## SRR10912089 Serous_EOC Control 3C
## SRR10912090 Serous_EOC Control 3C
## SRR10912091 Serous_EOC Control 3C
## SRR10912092 Serous_EOC Control 3C
## SRR10912093 Serous_EOC Control 4
## SRR10912094 Serous_EOC Control 3C
## SRR10912095 Serous_EOC Control 3C
## SRR10912096 Serous_EOC Control 3C
## SRR10912097 Serous_EOC Control 3C
## SRR10912098 Serous_EOC Control 3C
## SRR10912099 Serous_EOC Control 3C
## SRR10912100 Serous_EOC Control 3C
## SRR10912101 Serous_EOC Control 3C
## SRR10912102 Serous_EOC Control 3C
## SRR10912103 Serous_EOC Control 2C
## SRR10912104 Serous_EOC Control 3C
## SRR10912105 Serous_EOC Control 4
## SRR10912106 Serous_EOC Control 3C
## SRR10912107 Serous_EOC Control 3C
## SRR10912108 Serous_EOC Control 3C
## SRR10912109 Serous_EOC Control 4
## SRR10912110 Serous_EOC Control 3C
## SRR10912111 Serous_EOC Control 3C
## SRR10912112 Serous_EOC Control 3C
## SRR10912113 Serous_EOC Control 3C
## SRR10912114 Serous_EOC Control 3C
## SRR10912115 Serous_EOC Control 4
## SRR10912116 Serous_EOC Control 3C
## SRR10912117 Serous_EOC Control 3C
## SRR10912118 Serous_EOC Control 3C
## SRR10912119 Serous_EOC Control 3C
## SRR10912120 Serous_EOC Control 4
## SRR10912121 Serous_EOC Control 3C
## SRR10912122 Serous_EOC Control 3C
## SRR10912123 Serous_EOC Control 3C
## SRR10912124 Serous_EOC Control 4
## SRR10912125 Serous_EOC Control 3C
## SRR10912126 Serous_EOC Control 4
## SRR10912127 Benign_Tissue Control patientid: 61
## SRR10912128 Benign_Tissue Control patientid: 62
## SRR10912129 Benign_Tissue Control patientid: 63
## SRR10912130 Benign_Tissue Control patientid: 64
## SRR10912131 Benign_Tissue Control patientid: 65
## SRR10912132 Benign_Tissue Control patientid: 66
## SRR10912133 Benign_Tissue Control patientid: 67
## SRR10912134 Benign_Tissue Control patientid: 68
## SRR10912135 Benign_Tissue Control patientid: 69
## SRR10912136 Benign_Tissue Control patientid: 70
## SRR10912137 Benign_Tissue Control patientid: 71
## SRR10912138 Ascites Control 4A
## SRR10912139 Ascites Control 3C
## SRR10912140 Ascites Control 3C
## SRR10912141 Ascites Control 3C
## SRR10912142 Serous_EOC Control 4A
## SRR10912143 Serous_EOC Control 3C
## SRR10912144 Benign_Tissue Control patientid: 78
## SRR10912145 Serous_EOC Control 3C
## SRR10912146 Serous_EOC Control 3C
## SRR10912147 Ascites Control 3C
## SRR10912148 Ascites Control 3C
## SRR10912149 Ascites Control 3C
## SRR10912150 Ascites Control 3C
## SRR10912151 Ascites Control 3C
## SRR10912152 Ascites Control 4
## SRR10912153 Ascites Control 3C
## SRR10912154 Ascites Control 3C
## SRR10912155 Ascites Control 3C
## SRR10912156 Ascites Control 3C
## SRR10912157 Ascites Control 3C
## SRR10912158 Ascites Control 3C
## SRR10912159 Ascites Control 3C
## SRR10912160 Ascites Control 3C
## SRR10912161 Ascites Control 2C
## SRR10912162 Ascites Control 3C
## SRR10912163 Ascites Control 4
## SRR10912164 Ascites Control 3C
## SRR10912165 Ascites Control 3C
## SRR10912166 Ascites Control 3C
## SRR10912167 Ascites Control 3C
## SRR10912168 Ascites Control 4
## SRR10912169 Ascites Control 3C
## SRR10912170 Ascites Control 3C
## SRR10912171 Ascites Control 3C
## SRR10912172 Ascites Control 4This takes the metadata, or column data, file and corrects the columns and index of the dataframe, then converts the file from a data frame to a matrix.
Check if rownames and columnnames Match
all(rownames(column_data) %in% colnames(cnts))## [1] TRUEall(rownames(column_data) == colnames(cnts))## [1] TRUEThis checks if the rownames and column names match.
DeSeq2
library("DESeq2")## Loading required package: S4Vectors## Loading required package: stats4## Loading required package: BiocGenerics##
## Attaching package: 'BiocGenerics'## The following objects are masked from 'package:stats':
##
## IQR, mad, sd, var, xtabs## The following objects are masked from 'package:base':
##
## anyDuplicated, aperm, append, as.data.frame, basename, cbind,
## colnames, dirname, do.call, duplicated, eval, evalq, Filter, Find,
## get, grep, grepl, intersect, is.unsorted, lapply, Map, mapply,
## match, mget, order, paste, pmax, pmax.int, pmin, pmin.int,
## Position, rank, rbind, Reduce, rownames, sapply, setdiff, sort,
## table, tapply, union, unique, unsplit, which.max, which.min##
## Attaching package: 'S4Vectors'## The following object is masked from 'package:utils':
##
## findMatches## The following objects are masked from 'package:base':
##
## expand.grid, I, unname## Loading required package: IRanges## Loading required package: GenomicRanges## Loading required package: GenomeInfoDb## Loading required package: SummarizedExperiment## Loading required package: MatrixGenerics## Loading required package: matrixStats##
## Attaching package: 'MatrixGenerics'## The following objects are masked from 'package:matrixStats':
##
## colAlls, colAnyNAs, colAnys, colAvgsPerRowSet, colCollapse,
## colCounts, colCummaxs, colCummins, colCumprods, colCumsums,
## colDiffs, colIQRDiffs, colIQRs, colLogSumExps, colMadDiffs,
## colMads, colMaxs, colMeans2, colMedians, colMins, colOrderStats,
## colProds, colQuantiles, colRanges, colRanks, colSdDiffs, colSds,
## colSums2, colTabulates, colVarDiffs, colVars, colWeightedMads,
## colWeightedMeans, colWeightedMedians, colWeightedSds,
## colWeightedVars, rowAlls, rowAnyNAs, rowAnys, rowAvgsPerColSet,
## rowCollapse, rowCounts, rowCummaxs, rowCummins, rowCumprods,
## rowCumsums, rowDiffs, rowIQRDiffs, rowIQRs, rowLogSumExps,
## rowMadDiffs, rowMads, rowMaxs, rowMeans2, rowMedians, rowMins,
## rowOrderStats, rowProds, rowQuantiles, rowRanges, rowRanks,
## rowSdDiffs, rowSds, rowSums2, rowTabulates, rowVarDiffs, rowVars,
## rowWeightedMads, rowWeightedMeans, rowWeightedMedians,
## rowWeightedSds, rowWeightedVars## Loading required package: Biobase## Welcome to Bioconductor
##
## Vignettes contain introductory material; view with
## 'browseVignettes()'. To cite Bioconductor, see
## 'citation("Biobase")', and for packages 'citation("pkgname")'.##
## Attaching package: 'Biobase'## The following object is masked from 'package:MatrixGenerics':
##
## rowMedians## The following objects are masked from 'package:matrixStats':
##
## anyMissing, rowMediansdds_Tissue <- DESeqDataSetFromMatrix(countData = cnts,
colData = column_data,
design = ~ Tissue)## Warning in DESeqDataSet(se, design = design, ignoreRank): some variables in
## design formula are characters, converting to factorsdds_Tissue## class: DESeqDataSet
## dim: 20643 121
## metadata(1): version
## assays(1): counts
## rownames(20643): ENSG00000000003 ENSG00000000005 ... ENSG00000273439
## ENSG00000273452
## rowData names(0):
## colnames(121): SRR10912052 SRR10912053 ... SRR10912171 SRR10912172
## colData names(3): Tissue Treatment Stagedds_Treatment <- DESeqDataSetFromMatrix(countData = cnts,
colData = column_data,
design = ~ Treatment)
dds_Treatment## class: DESeqDataSet
## dim: 20643 121
## metadata(1): version
## assays(1): counts
## rownames(20643): ENSG00000000003 ENSG00000000005 ... ENSG00000273439
## ENSG00000273452
## rowData names(0):
## colnames(121): SRR10912052 SRR10912053 ... SRR10912171 SRR10912172
## colData names(3): Tissue Treatment StagePre-filtering
keep <- rowSums(counts(dds_Tissue)) >= 10
dds_Tissue <- dds_Tissue[keep,]
keep <- rowSums(counts(dds_Treatment)) >= 10
dds_Treatment <- dds_Treatment[keep,]Prefiltering can improve visualization as features with no information are not plotted.
Factor Levels
dds_Tissue$Tissue <- relevel(dds_Tissue$Tissue, ref = "Benign_Tissue")
dds_Treatment$Treatment <- factor(dds_Treatment$Treatment, levels = c("Control","Platinum"))“By default, R will choose a reference level for factors
based on alphabetical order. Then, if you never tell the DESeq2
functions which level you want to compare against (e.g. which level
represents the control group), the comparisons will be based on the
alphabetical order of the levels. There are two solutions: you can
either explicitly tell results which comparison to make using
the contrast argument (this will be shown later), or you
can explicitly set the factors levels. In order to see the change of
reference levels reflected in the results names, you need to either run
DESeq or nbinomWaldTest/nbinomLRT
after the re-leveling operation.”
Setting Benign Tissue as the reference level.
Differential expression analysis
dds_Tissue <- DESeq(dds_Tissue)## estimating size factors## estimating dispersions## gene-wise dispersion estimates## mean-dispersion relationship## final dispersion estimates## fitting model and testing## -- replacing outliers and refitting for 1168 genes
## -- DESeq argument 'minReplicatesForReplace' = 7
## -- original counts are preserved in counts(dds)## estimating dispersions## fitting model and testingresultsNames(dds_Tissue)## [1] "Intercept" "Tissue_Ascites_vs_Benign_Tissue"
## [3] "Tissue_Serous_EOC_vs_Benign_Tissue"dds_Treatment <- DESeq(dds_Treatment)## estimating size factors## estimating dispersions## gene-wise dispersion estimates## mean-dispersion relationship## final dispersion estimates## fitting model and testing## -- replacing outliers and refitting for 1583 genes
## -- DESeq argument 'minReplicatesForReplace' = 7
## -- original counts are preserved in counts(dds)## estimating dispersions## fitting model and testingresultsNames(dds_Treatment)## [1] "Intercept" "Treatment_Platinum_vs_Control"#summary(res)res_Serous_EOC <- results(dds_Tissue, name="Tissue_Serous_EOC_vs_Benign_Tissue")
res_Serous_EOC## log2 fold change (MLE): Tissue Serous EOC vs Benign Tissue
## Wald test p-value: Tissue Serous EOC vs Benign Tissue
## DataFrame with 19975 rows and 6 columns
## baseMean log2FoldChange lfcSE stat pvalue
## <numeric> <numeric> <numeric> <numeric> <numeric>
## ENSG00000000003 1359.4263 -0.790170 0.248072 -3.185252 1.44628e-03
## ENSG00000000005 28.6264 -0.512972 0.913941 -0.561275 5.74610e-01
## ENSG00000000419 945.1298 0.408815 0.156660 2.609561 9.06584e-03
## ENSG00000000457 339.0323 -0.139602 0.124784 -1.118752 2.63246e-01
## ENSG00000000460 189.3951 1.346737 0.231249 5.823762 5.75375e-09
## ... ... ... ... ... ...
## ENSG00000273294 18.6100019 0.115738 0.774778 0.1493822 0.881252
## ENSG00000273331 1.6135503 2.457218 1.142924 2.1499398 0.031560
## ENSG00000273398 5.7510137 0.562653 0.536606 1.0485403 0.294390
## ENSG00000273439 22.4666648 0.669630 0.415259 1.6125583 0.106841
## ENSG00000273452 0.0657423 0.197680 4.954781 0.0398967 0.968175
## padj
## <numeric>
## ENSG00000000003 3.91623e-03
## ENSG00000000005 6.64751e-01
## ENSG00000000419 1.99003e-02
## ENSG00000000457 3.58693e-01
## ENSG00000000460 5.14674e-08
## ... ...
## ENSG00000273294 0.9141591
## ENSG00000273331 0.0593248
## ENSG00000273398 0.3925197
## ENSG00000273439 0.1682453
## ENSG00000273452 NAresLFC_Serous_EOC <- lfcShrink(dds_Tissue, coef="Tissue_Serous_EOC_vs_Benign_Tissue", type="apeglm")## using 'apeglm' for LFC shrinkage. If used in published research, please cite:
## Zhu, A., Ibrahim, J.G., Love, M.I. (2018) Heavy-tailed prior distributions for
## sequence count data: removing the noise and preserving large differences.
## Bioinformatics. https://doi.org/10.1093/bioinformatics/bty895resLFC_Serous_EOC## log2 fold change (MAP): Tissue Serous EOC vs Benign Tissue
## Wald test p-value: Tissue Serous EOC vs Benign Tissue
## DataFrame with 19975 rows and 5 columns
## baseMean log2FoldChange lfcSE pvalue padj
## <numeric> <numeric> <numeric> <numeric> <numeric>
## ENSG00000000003 1359.4263 -0.707614 0.238990 1.44628e-03 3.91623e-03
## ENSG00000000005 28.6264 -0.256894 0.633896 5.74610e-01 6.64751e-01
## ENSG00000000419 945.1298 0.397866 0.155273 9.06584e-03 1.99003e-02
## ENSG00000000457 339.0323 -0.136250 0.123498 2.63246e-01 3.58693e-01
## ENSG00000000460 189.3951 1.301129 0.234105 5.75375e-09 5.14674e-08
## ... ... ... ... ... ...
## ENSG00000273294 18.6100019 0.0682127 0.587723 0.881252 0.9141591
## ENSG00000273331 1.6135503 3.6905904 1.538611 0.031560 0.0593248
## ENSG00000273398 5.7510137 0.4212642 0.487248 0.294390 0.3925197
## ENSG00000273439 22.4666648 0.5644270 0.401099 0.106841 0.1682453
## ENSG00000273452 0.0657423 0.0337805 0.882672 0.968175 NAres_Ascites <- results(dds_Tissue, name="Tissue_Ascites_vs_Benign_Tissue")
res_Ascites## log2 fold change (MLE): Tissue Ascites vs Benign Tissue
## Wald test p-value: Tissue Ascites vs Benign Tissue
## DataFrame with 19975 rows and 6 columns
## baseMean log2FoldChange lfcSE stat pvalue
## <numeric> <numeric> <numeric> <numeric> <numeric>
## ENSG00000000003 1359.4263 -0.914937 0.275916 -3.31600 9.13167e-04
## ENSG00000000005 28.6264 -4.565551 1.029674 -4.43398 9.25110e-06
## ENSG00000000419 945.1298 0.487980 0.175152 2.78604 5.33562e-03
## ENSG00000000457 339.0323 -0.448853 0.141450 -3.17321 1.50761e-03
## ENSG00000000460 189.3951 1.099443 0.259719 4.23321 2.30384e-05
## ... ... ... ... ... ...
## ENSG00000273294 18.6100019 -0.553297 0.866688 -0.638404 0.5232107
## ENSG00000273331 1.6135503 2.834898 1.244909 2.277193 0.0227747
## ENSG00000273398 5.7510137 -0.419677 0.607245 -0.691117 0.4894919
## ENSG00000273439 22.4666648 -0.072854 0.470010 -0.155005 0.8768173
## ENSG00000273452 0.0657423 2.226346 5.458155 0.407893 0.6833519
## padj
## <numeric>
## ENSG00000000003 2.17650e-03
## ENSG00000000005 3.16468e-05
## ENSG00000000419 1.08888e-02
## ENSG00000000457 3.44937e-03
## ENSG00000000460 7.32642e-05
## ... ...
## ENSG00000273294 0.6021114
## ENSG00000273331 0.0401088
## ENSG00000273398 0.5705074
## ENSG00000273439 0.9009092
## ENSG00000273452 0.7426654resLFC_Ascites <- lfcShrink(dds_Tissue, coef="Tissue_Ascites_vs_Benign_Tissue", type="apeglm")## using 'apeglm' for LFC shrinkage. If used in published research, please cite:
## Zhu, A., Ibrahim, J.G., Love, M.I. (2018) Heavy-tailed prior distributions for
## sequence count data: removing the noise and preserving large differences.
## Bioinformatics. https://doi.org/10.1093/bioinformatics/bty895resLFC_Ascites## log2 fold change (MAP): Tissue Ascites vs Benign Tissue
## Wald test p-value: Tissue Ascites vs Benign Tissue
## DataFrame with 19975 rows and 5 columns
## baseMean log2FoldChange lfcSE pvalue padj
## <numeric> <numeric> <numeric> <numeric> <numeric>
## ENSG00000000003 1359.4263 -0.837964 0.267011 9.13167e-04 2.17650e-03
## ENSG00000000005 28.6264 -0.844550 0.833044 9.25110e-06 3.16468e-05
## ENSG00000000419 945.1298 0.475800 0.173115 5.33562e-03 1.08888e-02
## ENSG00000000457 339.0323 -0.441299 0.139981 1.50761e-03 3.44937e-03
## ENSG00000000460 189.3951 1.050195 0.259593 2.30384e-05 7.32642e-05
## ... ... ... ... ... ...
## ENSG00000273294 18.6100019 -0.4558388 0.684649 0.5232107 0.6021114
## ENSG00000273331 1.6135503 0.8943902 1.171554 0.0227747 0.0401088
## ENSG00000273398 5.7510137 -0.5741662 0.558885 0.4894919 0.5705074
## ENSG00000273439 22.4666648 -0.1303708 0.426959 0.8768173 0.9009092
## ENSG00000273452 0.0657423 0.0804353 1.022513 0.6833519 0.7426654res_Treatment <- results(dds_Treatment, name="Treatment_Platinum_vs_Control")
res_Treatment## log2 fold change (MLE): Treatment Platinum vs Control
## Wald test p-value: Treatment Platinum vs Control
## DataFrame with 19975 rows and 6 columns
## baseMean log2FoldChange lfcSE stat pvalue
## <numeric> <numeric> <numeric> <numeric> <numeric>
## ENSG00000000003 1359.4263 -0.3085477 0.213204 -1.447194 0.147843
## ENSG00000000005 32.3854 0.3109034 0.842892 0.368853 0.712237
## ENSG00000000419 945.1298 0.0916123 0.132497 0.691430 0.489295
## ENSG00000000457 339.0323 0.0517086 0.107236 0.482192 0.629669
## ENSG00000000460 189.3951 -0.2795861 0.207999 -1.344168 0.178894
## ... ... ... ... ... ...
## ENSG00000273294 18.6100019 -0.935388 0.638808 -1.464269 0.1431203
## ENSG00000273331 1.6135503 0.387927 0.858770 0.451724 0.6514677
## ENSG00000273398 5.7510137 -0.334968 0.456960 -0.733036 0.4635363
## ENSG00000273439 22.4666648 -0.737267 0.344510 -2.140043 0.0323513
## ENSG00000273452 0.0657423 -0.698464 4.087144 -0.170893 0.8643080
## padj
## <numeric>
## ENSG00000000003 0.294437
## ENSG00000000005 0.831205
## ENSG00000000419 0.657758
## ENSG00000000457 0.771404
## ENSG00000000460 0.336298
## ... ...
## ENSG00000273294 0.2878137
## ENSG00000273331 0.7883023
## ENSG00000273398 0.6364602
## ENSG00000273439 0.0978889
## ENSG00000273452 NAresLFC_Treatment <- lfcShrink(dds_Treatment, coef="Treatment_Platinum_vs_Control", type="apeglm")## using 'apeglm' for LFC shrinkage. If used in published research, please cite:
## Zhu, A., Ibrahim, J.G., Love, M.I. (2018) Heavy-tailed prior distributions for
## sequence count data: removing the noise and preserving large differences.
## Bioinformatics. https://doi.org/10.1093/bioinformatics/bty895resLFC_Treatment## log2 fold change (MAP): Treatment Platinum vs Control
## Wald test p-value: Treatment Platinum vs Control
## DataFrame with 19975 rows and 5 columns
## baseMean log2FoldChange lfcSE pvalue padj
## <numeric> <numeric> <numeric> <numeric> <numeric>
## ENSG00000000003 1359.4263 -6.55471e-06 0.00144267 0.147843 0.294437
## ENSG00000000005 32.3854 5.52645e-07 0.00144269 0.712237 0.831205
## ENSG00000000419 945.1298 4.03178e-05 0.00144289 0.489295 0.657758
## ENSG00000000457 339.0323 4.74106e-06 0.00144257 0.629669 0.771404
## ENSG00000000460 189.3951 -6.27916e-06 0.00144267 0.178894 0.336298
## ... ... ... ... ... ...
## ENSG00000273294 18.6100019 -1.85032e-06 0.00144269 0.1431203 0.2878137
## ENSG00000273331 1.6135503 3.60722e+00 0.81074595 0.6514677 0.7883023
## ENSG00000273398 5.7510137 -1.32970e-06 0.00144269 0.4635363 0.6364602
## ENSG00000273439 22.4666648 -5.40768e-06 0.00144269 0.0323513 0.0978889
## ENSG00000273452 0.0657423 -1.05593e-07 0.00144270 0.8643080 NAUsing lfcShrik to make vizualization and ranking of genes better.
resLFCOrdered_Serous_EOC <- resLFC_Serous_EOC[order(resLFC_Serous_EOC$pvalue),]
resLFCOrdered_Serous_EOC## log2 fold change (MAP): Tissue Serous EOC vs Benign Tissue
## Wald test p-value: Tissue Serous EOC vs Benign Tissue
## DataFrame with 19975 rows and 5 columns
## baseMean log2FoldChange lfcSE pvalue padj
## <numeric> <numeric> <numeric> <numeric> <numeric>
## ENSG00000100380 5724.9650 -2.36653 0.160123 6.74717e-50 9.08323e-46
## ENSG00000154545 61.4084 11.23218 3.529445 9.46367e-50 9.08323e-46
## ENSG00000100227 1539.3231 -1.92776 0.143622 2.46081e-41 1.57459e-37
## ENSG00000139734 173.2278 4.27102 0.321202 9.45548e-41 4.53768e-37
## ENSG00000069966 566.3977 -2.15348 0.163485 4.28627e-40 1.64559e-36
## ... ... ... ... ... ...
## ENSG00000184007 3630.87 0.000342523 0.127981 0.999963 0.999963
## ENSG00000154537 0.00 0.069290955 0.887184 1.000000 NA
## ENSG00000242366 0.00 0.061543666 0.885882 1.000000 NA
## ENSG00000255863 0.00 0.071625991 0.887520 1.000000 NA
## ENSG00000268485 0.00 0.068737657 0.887100 1.000000 NAresLFCOrdered_Ascites <- resLFC_Ascites[order(resLFC_Ascites$pvalue),]
resLFCOrdered_Ascites## log2 fold change (MAP): Tissue Ascites vs Benign Tissue
## Wald test p-value: Tissue Ascites vs Benign Tissue
## DataFrame with 19975 rows and 5 columns
## baseMean log2FoldChange lfcSE pvalue padj
## <numeric> <numeric> <numeric> <numeric> <numeric>
## ENSG00000115461 16817.080 -6.76674 0.420994 4.15654e-58 8.30102e-54
## ENSG00000100380 5724.965 -2.74955 0.178141 2.86234e-54 2.85819e-50
## ENSG00000072840 769.001 -3.36303 0.223246 7.75342e-52 5.16145e-48
## ENSG00000100227 1539.323 -2.41606 0.161322 2.93616e-51 1.46595e-47
## ENSG00000113658 1041.526 -2.28680 0.158517 1.05393e-47 4.20962e-44
## ... ... ... ... ... ...
## ENSG00000147596 3.33851 -0.1884468 0.527538 0.999197 0.999197
## ENSG00000154537 0.00000 -0.0302969 1.017772 1.000000 NA
## ENSG00000242366 0.00000 -0.0277997 1.017821 1.000000 NA
## ENSG00000255863 0.00000 -0.0311519 1.017759 1.000000 NA
## ENSG00000268485 0.00000 -0.0301904 1.017780 1.000000 NAresLFC_Treatment <- resLFC_Treatment[order(resLFC_Treatment$pvalue),]
resLFC_Treatment## log2 fold change (MAP): Treatment Platinum vs Control
## Wald test p-value: Treatment Platinum vs Control
## DataFrame with 19975 rows and 5 columns
## baseMean log2FoldChange lfcSE pvalue padj
## <numeric> <numeric> <numeric> <numeric> <numeric>
## ENSG00000232196 2250.5603 -9.90611e-06 0.00144271 2.95619e-39 5.78734e-35
## ENSG00000152954 427.3347 -6.07833e+00 0.48871479 6.68560e-38 6.54420e-34
## ENSG00000214954 473.4323 -5.54753e+00 0.45616227 1.91708e-28 1.25102e-24
## ENSG00000185873 27.7571 -9.50079e-06 0.00144271 9.82233e-28 4.80730e-24
## ENSG00000130032 1514.0007 -9.76177e-06 0.00144271 2.23886e-27 8.76602e-24
## ... ... ... ... ... ...
## ENSG00000255863 0 -1.06574e-07 0.0014427 1 NA
## ENSG00000258389 0 5.41012e-07 0.0014427 1 NA
## ENSG00000261833 0 -1.02698e-07 0.0014427 1 NA
## ENSG00000268427 0 -1.02498e-07 0.0014427 1 NA
## ENSG00000268485 0 5.98608e-07 0.0014427 1 NAOrdered by p-value.
summary(resLFCOrdered_Serous_EOC)##
## out of 19971 with nonzero total read count
## adjusted p-value < 0.1
## LFC > 0 (up) : 6618, 33%
## LFC < 0 (down) : 4518, 23%
## outliers [1] : 0, 0%
## low counts [2] : 779, 3.9%
## (mean count < 0)
## [1] see 'cooksCutoff' argument of ?results
## [2] see 'independentFiltering' argument of ?resultssummary(resLFCOrdered_Ascites)##
## out of 19971 with nonzero total read count
## adjusted p-value < 0.1
## LFC > 0 (up) : 6636, 33%
## LFC < 0 (down) : 6123, 31%
## outliers [1] : 0, 0%
## low counts [2] : 4, 0.02%
## (mean count < 0)
## [1] see 'cooksCutoff' argument of ?results
## [2] see 'independentFiltering' argument of ?resultssummary(resLFC_Treatment)##
## out of 19965 with nonzero total read count
## adjusted p-value < 0.1
## LFC > 0 (up) : 2520, 13%
## LFC < 0 (down) : 3997, 20%
## outliers [1] : 0, 0%
## low counts [2] : 398, 2%
## (mean count < 0)
## [1] see 'cooksCutoff' argument of ?results
## [2] see 'independentFiltering' argument of ?resultsMA Plots
plotMA(resLFCOrdered_Serous_EOC, ylim=c(-2,2))
idx_LFCOrdered_Serous_EOC <- identify(resLFCOrdered_Serous_EOC$baseMean, resLFCOrdered_Serous_EOC$log2FoldChange)
rownames(resLFCOrdered_Serous_EOC)[idx_LFCOrdered_Serous_EOC]## character(0)plotMA(resLFCOrdered_Ascites, ylim=c(-2,2))
plotMA(resLFC_Treatment, ylim=c(-2,2))
plotMA(resLFCOrdered_Serous_EOC, ylim=c(-2,2))
idx_LFCOrdered_Serous_EOC <- identify(resLFCOrdered_Serous_EOC$baseMean, resLFCOrdered_Serous_EOC$log2FoldChange)
rownames(resLFCOrdered_Serous_EOC)[idx_LFCOrdered_Serous_EOC]## character(0)resNorm <- lfcShrink(dds_Tissue, coef=3, type="normal")## using 'normal' for LFC shrinkage, the Normal prior from Love et al (2014).
##
## Note that type='apeglm' and type='ashr' have shown to have less bias than type='normal'.
## See ?lfcShrink for more details on shrinkage type, and the DESeq2 vignette.
## Reference: https://doi.org/10.1093/bioinformatics/bty895resAsh <- lfcShrink(dds_Tissue, coef=3, type="ashr")## using 'ashr' for LFC shrinkage. If used in published research, please cite:
## Stephens, M. (2016) False discovery rates: a new deal. Biostatistics, 18:2.
## https://doi.org/10.1093/biostatistics/kxw041resLFC <- lfcShrink(dds_Tissue, coef=3, type="apeglm")## using 'apeglm' for LFC shrinkage. If used in published research, please cite:
## Zhu, A., Ibrahim, J.G., Love, M.I. (2018) Heavy-tailed prior distributions for
## sequence count data: removing the noise and preserving large differences.
## Bioinformatics. https://doi.org/10.1093/bioinformatics/bty895par(mfrow=c(1,3), mar=c(4,4,2,1))
xlim <- c(1,1e5); ylim <- c(-3,3)
plotMA(resLFC, xlim=xlim, ylim=ylim, main="apeglm")
plotMA(resNorm, xlim=xlim, ylim=ylim, main="normal")
plotMA(resAsh, xlim=xlim, ylim=ylim, main="ashr")
resNorm <- lfcShrink(dds_Treatment, coef=2, type="normal")## using 'normal' for LFC shrinkage, the Normal prior from Love et al (2014).
##
## Note that type='apeglm' and type='ashr' have shown to have less bias than type='normal'.
## See ?lfcShrink for more details on shrinkage type, and the DESeq2 vignette.
## Reference: https://doi.org/10.1093/bioinformatics/bty895resAsh <- lfcShrink(dds_Treatment, coef=2, type="ashr")## using 'ashr' for LFC shrinkage. If used in published research, please cite:
## Stephens, M. (2016) False discovery rates: a new deal. Biostatistics, 18:2.
## https://doi.org/10.1093/biostatistics/kxw041resLFC <- lfcShrink(dds_Treatment, coef=2, type="apeglm")## using 'apeglm' for LFC shrinkage. If used in published research, please cite:
## Zhu, A., Ibrahim, J.G., Love, M.I. (2018) Heavy-tailed prior distributions for
## sequence count data: removing the noise and preserving large differences.
## Bioinformatics. https://doi.org/10.1093/bioinformatics/bty895par(mfrow=c(1,3), mar=c(4,4,2,1))
xlim <- c(1,1e5); ylim <- c(-3,3)
plotMA(resLFC, xlim=xlim, ylim=ylim, main="apeglm")
plotMA(resNorm, xlim=xlim, ylim=ylim, main="normal")
plotMA(resAsh, xlim=xlim, ylim=ylim, main="ashr")
Plot Counts
plotCounts(dds_Tissue, gene=which.min(resLFC_Serous_EOC$padj), intgroup="Tissue")
plotCounts(dds_Tissue, gene=which.min(resLFC_Ascites$padj), intgroup="Tissue")
plotCounts(dds_Treatment, gene=which.min(resLFC_Treatment$padj), intgroup="Treatment")
plotCounts(dds_Tissue, gene=which.max(resLFC_Serous_EOC$padj), intgroup="Tissue")
plotCounts(dds_Tissue, gene=which.max(resLFC_Ascites$padj), intgroup="Tissue")
plotCounts(dds_Treatment, gene=which.max(resLFC_Treatment$padj), intgroup="Treatment")
Data Transformations + Variance
vsd <- vst(dds_Tissue, blind=FALSE)
#rld <- rlog(dds_Tissue, blind=FALSE)
ntd <- normTransform(dds_Tissue)
library("vsn")
head(assay(vsd), 10)## SRR10912052 SRR10912053 SRR10912054 SRR10912055 SRR10912056
## ENSG00000000003 9.681800 11.409086 10.373103 10.861253 9.804261
## ENSG00000000005 5.657905 7.242788 7.324413 5.657905 5.657905
## ENSG00000000419 10.143497 9.959867 10.159948 10.482950 10.201271
## ENSG00000000457 8.778506 8.722905 8.701035 8.582886 8.955522
## ENSG00000000460 8.362459 7.709567 8.302766 7.852794 7.959261
## ENSG00000000938 8.780148 8.866326 8.220189 6.892303 8.679407
## ENSG00000000971 11.347021 12.181838 12.121527 9.671775 13.727215
## ENSG00000001036 10.324587 10.764334 10.538936 9.731350 10.829804
## ENSG00000001084 8.812592 9.807770 9.398510 9.429871 9.110757
## ENSG00000001167 9.789864 9.155460 9.201835 9.320337 9.062793
## SRR10912057 SRR10912058 SRR10912059 SRR10912060 SRR10912061
## ENSG00000000003 10.682838 9.112636 8.461819 10.639874 10.202538
## ENSG00000000005 5.834671 6.058665 5.657905 6.745606 5.657905
## ENSG00000000419 10.023898 9.791563 9.282585 10.369044 9.815419
## ENSG00000000457 8.986655 8.742551 8.358014 8.744751 8.705403
## ENSG00000000460 8.710638 7.591303 7.993237 7.186378 5.657905
## ENSG00000000938 7.648104 10.074286 10.837538 7.853983 8.663492
## ENSG00000000971 11.194760 13.162126 11.148968 12.214480 11.413431
## ENSG00000001036 10.445657 10.926095 11.153241 10.424914 10.027705
## ENSG00000001084 9.669354 9.722022 9.641572 8.383807 8.633388
## ENSG00000001167 9.844475 9.114695 8.595770 9.295579 9.002105
## SRR10912062 SRR10912063 SRR10912064 SRR10912065 SRR10912066
## ENSG00000000003 10.479004 9.696886 11.402766 10.635603 10.662796
## ENSG00000000005 5.657905 6.540222 5.657905 6.048287 6.043689
## ENSG00000000419 9.998549 9.591109 9.994770 9.665286 10.498476
## ENSG00000000457 7.540698 8.059312 9.196570 9.322005 8.367355
## ENSG00000000460 8.117656 8.041946 8.708845 8.244808 8.577920
## ENSG00000000938 8.856192 9.497592 8.243940 8.463016 6.875277
## ENSG00000000971 12.171610 12.334648 12.146060 12.311195 10.779275
## ENSG00000001036 10.932153 10.660764 10.754683 10.506151 10.429774
## ENSG00000001084 9.437551 10.370967 9.628905 9.583375 9.188028
## ENSG00000001167 9.064419 9.085622 9.245514 9.145594 9.537424
## SRR10912067 SRR10912068 SRR10912069 SRR10912070 SRR10912071
## ENSG00000000003 9.514325 9.748332 8.928942 10.655050 9.182824
## ENSG00000000005 6.062590 6.024113 6.033176 5.906723 5.966233
## ENSG00000000419 10.218577 10.090740 9.729964 10.343107 10.317958
## ENSG00000000457 8.856011 8.879201 8.956091 8.810905 8.523476
## ENSG00000000460 7.857643 8.024855 7.788848 8.499530 7.790599
## ENSG00000000938 11.264581 9.000551 10.138362 7.652670 9.462310
## ENSG00000000971 10.352780 11.321204 11.565934 11.378536 13.652808
## ENSG00000001036 10.493664 10.356254 10.325404 10.304064 10.271805
## ENSG00000001084 10.077249 9.041806 9.421800 9.444026 9.161434
## ENSG00000001167 9.114372 9.555633 9.289950 9.149040 8.936867
## SRR10912072 SRR10912073 SRR10912074 SRR10912075 SRR10912076
## ENSG00000000003 9.789805 11.805167 10.369948 9.528750 9.369386
## ENSG00000000005 5.859382 6.659314 7.349422 5.657905 7.217420
## ENSG00000000419 9.869783 10.405773 10.355734 9.860681 9.958212
## ENSG00000000457 8.597413 9.256500 9.136575 9.553261 9.034634
## ENSG00000000460 8.500700 8.517670 8.866012 8.830151 7.911606
## ENSG00000000938 8.691349 6.881988 8.925536 8.722829 9.153147
## ENSG00000000971 12.707167 8.275037 10.985326 10.365997 12.363049
## ENSG00000001036 10.318118 9.742067 10.369948 9.985780 9.956461
## ENSG00000001084 9.152122 9.668891 9.778350 9.639703 9.747536
## ENSG00000001167 9.206526 10.435005 9.878489 9.632057 9.041300
## SRR10912077 SRR10912078 SRR10912079 SRR10912080 SRR10912081
## ENSG00000000003 9.237389 11.686149 11.652404 10.949399 9.884800
## ENSG00000000005 6.335134 7.790812 6.513903 5.657905 5.853438
## ENSG00000000419 9.695735 10.505656 10.398156 11.690905 10.842806
## ENSG00000000457 8.858825 8.695915 8.726438 8.675206 8.720754
## ENSG00000000460 7.905158 8.464077 7.902284 8.496152 7.790266
## ENSG00000000938 9.729510 7.692708 7.216529 7.834771 9.635456
## ENSG00000000971 12.557434 11.256344 12.572694 11.313913 12.410532
## ENSG00000001036 10.501549 10.498353 10.290944 11.849828 11.387249
## ENSG00000001084 9.420927 9.000285 9.335599 10.360706 9.871949
## ENSG00000001167 9.121185 10.372169 9.683824 9.518110 8.876175
## SRR10912082 SRR10912083 SRR10912084 SRR10912085 SRR10912086
## ENSG00000000003 10.731144 10.626092 10.391571 9.856398 10.063216
## ENSG00000000005 5.657905 5.657905 5.657905 8.312581 9.890164
## ENSG00000000419 10.106190 10.459343 10.316807 10.277954 9.932346
## ENSG00000000457 9.471265 8.715264 8.799716 8.959380 8.697414
## ENSG00000000460 9.322198 8.184746 8.685387 7.888386 7.847349
## ENSG00000000938 8.230518 7.908309 8.739724 8.987737 9.060709
## ENSG00000000971 11.376767 12.608898 11.601021 13.279539 12.940010
## ENSG00000001036 9.516211 10.019321 10.552999 10.455881 10.626832
## ENSG00000001084 9.980664 9.376255 9.843706 9.739313 9.748638
## ENSG00000001167 10.132336 9.979692 9.969251 9.828365 8.866547
## SRR10912087 SRR10912088 SRR10912089 SRR10912090 SRR10912091
## ENSG00000000003 10.120827 10.632332 10.750688 10.123521 9.854531
## ENSG00000000005 8.319678 6.594620 6.053124 5.839106 6.584189
## ENSG00000000419 10.170484 10.062837 9.774579 9.751626 9.564163
## ENSG00000000457 9.112089 9.391028 8.996665 8.632326 8.715865
## ENSG00000000460 7.547950 8.538632 9.095365 7.636814 7.694398
## ENSG00000000938 10.262584 7.579326 7.497154 7.862851 7.712670
## ENSG00000000971 12.348515 10.132630 9.212140 11.335545 12.672983
## ENSG00000001036 10.334697 10.124089 11.204638 10.019598 10.092598
## ENSG00000001084 10.053189 9.559337 9.471102 8.982612 8.662969
## ENSG00000001167 9.244730 9.388645 9.431611 8.882670 8.772389
## SRR10912092 SRR10912093 SRR10912094 SRR10912095 SRR10912096
## ENSG00000000003 9.690714 9.773876 9.579622 10.102560 10.589301
## ENSG00000000005 8.596961 5.657905 5.888589 5.909040 5.966965
## ENSG00000000419 9.457148 10.056925 9.664049 10.200022 9.631575
## ENSG00000000457 8.509949 9.444754 8.808824 8.746192 8.825123
## ENSG00000000460 7.657982 8.724068 8.342910 7.990393 8.378910
## ENSG00000000938 9.221176 7.118030 7.965176 7.540845 7.778890
## ENSG00000000971 12.888248 9.997007 8.607472 11.344746 10.803650
## ENSG00000001036 10.539513 9.980753 9.948034 10.236310 9.962584
## ENSG00000001084 9.564194 8.540289 9.090626 9.241295 9.513569
## ENSG00000001167 9.260549 9.147596 8.697005 9.533242 9.597184
## SRR10912097 SRR10912098 SRR10912099 SRR10912100 SRR10912101
## ENSG00000000003 11.917693 10.727720 11.769118 10.730239 10.127653
## ENSG00000000005 6.775833 5.840411 5.657905 5.657905 5.657905
## ENSG00000000419 9.469495 9.354318 10.179571 9.974597 9.538607
## ENSG00000000457 8.889644 8.642477 9.182089 8.884587 8.232557
## ENSG00000000460 7.459067 7.679849 8.961568 8.881745 7.475655
## ENSG00000000938 6.978535 7.844592 6.463177 7.574433 7.505389
## ENSG00000000971 12.436809 10.611550 7.538526 10.600251 11.531618
## ENSG00000001036 11.046980 9.343513 10.749000 9.609189 9.802213
## ENSG00000001084 8.017013 8.703462 10.052257 9.155124 9.037496
## ENSG00000001167 9.456906 9.026969 10.028862 9.384264 9.795277
## SRR10912102 SRR10912103 SRR10912104 SRR10912105 SRR10912106
## ENSG00000000003 9.388029 10.768546 10.662914 10.441764 10.054127
## ENSG00000000005 8.013932 5.897503 6.325378 5.913889 6.284343
## ENSG00000000419 9.048870 11.036591 10.505685 9.626125 9.585915
## ENSG00000000457 8.615069 8.894731 8.733927 8.903725 8.604968
## ENSG00000000460 7.241883 9.250399 8.465158 7.981010 8.345257
## ENSG00000000938 9.351160 7.312097 8.292841 7.846434 7.805776
## ENSG00000000971 12.581507 8.971749 10.120269 9.127736 10.652037
## ENSG00000001036 10.216597 10.879143 10.023739 10.264039 9.860845
## ENSG00000001084 9.866794 9.628368 9.676119 10.030685 10.014120
## ENSG00000001167 8.979839 9.641085 9.898180 9.866292 9.464691
## SRR10912107 SRR10912108 SRR10912109 SRR10912110 SRR10912111
## ENSG00000000003 9.838507 9.999543 11.060003 10.957206 9.354553
## ENSG00000000005 6.397608 5.964393 6.371779 6.703043 5.977341
## ENSG00000000419 9.685932 9.556288 9.419346 9.590830 10.376476
## ENSG00000000457 8.182953 9.154727 8.625538 8.899580 8.966275
## ENSG00000000460 7.941163 8.197571 7.698895 8.043477 8.730311
## ENSG00000000938 8.566607 7.953876 7.563975 7.013731 7.760243
## ENSG00000000971 8.939289 10.077371 10.494718 9.782016 9.302175
## ENSG00000001036 9.663395 10.297584 10.618685 9.724768 10.291157
## ENSG00000001084 9.248075 9.322414 9.171858 10.042092 8.826827
## ENSG00000001167 8.959294 9.278970 9.378679 9.528836 9.958276
## SRR10912112 SRR10912113 SRR10912114 SRR10912115 SRR10912116
## ENSG00000000003 10.549125 10.169169 9.996266 11.065692 10.902300
## ENSG00000000005 5.657905 6.326860 6.035439 6.244687 5.657905
## ENSG00000000419 9.646125 10.161344 10.365855 10.019270 9.938859
## ENSG00000000457 9.053949 8.848088 8.596867 8.958203 8.876283
## ENSG00000000460 8.898865 8.496749 8.227468 8.459481 7.870837
## ENSG00000000938 7.434567 7.971937 7.364498 7.415554 7.100761
## ENSG00000000971 9.592380 10.305563 9.928723 10.170340 10.783769
## ENSG00000001036 10.057250 10.444526 10.124312 9.484405 9.825351
## ENSG00000001084 9.106949 9.324805 9.535727 8.553655 8.934686
## ENSG00000001167 9.230434 10.039492 10.144881 9.996473 9.352783
## SRR10912117 SRR10912118 SRR10912119 SRR10912120 SRR10912121
## ENSG00000000003 10.870113 9.791178 10.165048 10.569699 9.015377
## ENSG00000000005 5.857542 8.312728 6.100406 5.657905 5.657905
## ENSG00000000419 10.148289 8.973836 9.788153 10.042770 11.134821
## ENSG00000000457 8.950647 8.810180 8.608425 8.703128 8.566547
## ENSG00000000460 8.631807 7.705391 8.511551 8.292418 8.355014
## ENSG00000000938 7.420482 8.996380 8.176935 7.730380 8.085127
## ENSG00000000971 9.764560 11.878304 11.328730 11.961958 10.352330
## ENSG00000001036 10.424611 9.963537 10.097171 10.274575 10.198448
## ENSG00000001084 9.965713 9.548955 9.873345 8.995594 9.416324
## ENSG00000001167 9.134977 8.899180 10.474080 9.634069 8.879204
## SRR10912122 SRR10912123 SRR10912124 SRR10912125 SRR10912126
## ENSG00000000003 10.032690 10.087136 9.571070 10.769773 9.998517
## ENSG00000000005 5.917593 5.846894 5.657905 6.315303 6.083379
## ENSG00000000419 9.792948 10.325678 9.571070 9.705033 9.722388
## ENSG00000000457 8.428563 8.588788 8.249535 8.807468 8.599225
## ENSG00000000460 8.240850 8.653428 7.277281 8.350049 7.800779
## ENSG00000000938 7.679525 6.975433 10.694501 7.674585 8.072405
## ENSG00000000971 11.161452 11.393834 12.002600 11.354242 11.397260
## ENSG00000001036 10.541345 10.350856 10.290530 10.721072 10.140321
## ENSG00000001084 9.287902 9.326182 9.759159 9.412767 8.685759
## ENSG00000001167 9.587543 9.413703 9.824346 9.512946 8.729642
## SRR10912127 SRR10912128 SRR10912129 SRR10912130 SRR10912131
## ENSG00000000003 11.547610 11.424025 10.516266 11.203384 11.180243
## ENSG00000000005 7.042143 6.323793 7.599408 6.528323 6.551915
## ENSG00000000419 9.701540 9.541748 9.518293 9.632376 9.832567
## ENSG00000000457 9.114526 8.892260 8.684853 8.941324 8.812885
## ENSG00000000460 7.324757 7.774050 7.157379 7.359458 7.265680
## ENSG00000000938 6.771086 7.184867 7.623889 6.985543 6.601337
## ENSG00000000971 12.267779 12.696265 13.118564 12.085980 13.146398
## ENSG00000001036 10.692839 10.466388 10.706315 10.629241 10.610018
## ENSG00000001084 9.104741 9.686723 8.441543 8.613121 8.722025
## ENSG00000001167 9.628640 9.590581 9.300408 9.290722 9.664211
## SRR10912132 SRR10912133 SRR10912134 SRR10912135 SRR10912136
## ENSG00000000003 11.088528 11.352780 11.561319 11.162601 11.080903
## ENSG00000000005 8.238922 6.339365 6.550300 6.678743 7.238449
## ENSG00000000419 9.625747 9.802675 9.827910 9.769506 9.699860
## ENSG00000000457 8.813425 9.015513 9.136782 9.105001 8.784260
## ENSG00000000460 7.458393 7.507021 7.670911 7.478916 7.214206
## ENSG00000000938 8.016745 6.866548 6.405736 6.888255 6.576864
## ENSG00000000971 12.678047 13.099437 12.717134 12.577411 12.899147
## ENSG00000001036 10.349596 10.100810 10.584969 10.678041 10.562391
## ENSG00000001084 8.974784 8.976818 9.283510 8.705346 8.533921
## ENSG00000001167 9.426759 9.643702 9.801436 9.741259 9.475391
## SRR10912137 SRR10912138 SRR10912139 SRR10912140 SRR10912141
## ENSG00000000003 11.203735 10.462749 10.558232 9.352170 10.087283
## ENSG00000000005 7.240333 5.657905 5.876671 5.813607 5.927808
## ENSG00000000419 9.615981 10.624721 10.246869 9.749570 10.352192
## ENSG00000000457 8.866021 8.297502 8.579115 8.191324 8.416318
## ENSG00000000460 7.369684 7.890839 8.262644 7.675062 8.050090
## ENSG00000000938 6.577731 9.096736 8.967233 9.853455 9.511546
## ENSG00000000971 12.715642 11.472347 9.754289 11.011983 10.036134
## ENSG00000001036 10.475592 10.888224 10.426827 10.865948 9.660181
## ENSG00000001084 8.320703 8.781438 9.308742 9.567728 9.451031
## ENSG00000001167 9.782998 9.436212 9.320109 9.178360 9.010589
## SRR10912142 SRR10912143 SRR10912144 SRR10912145 SRR10912146
## ENSG00000000003 9.486470 11.422991 10.814979 9.913147 10.431310
## ENSG00000000005 6.957346 6.033677 5.657905 5.727410 6.085379
## ENSG00000000419 9.403874 10.087039 9.755601 9.566692 10.049713
## ENSG00000000457 8.619509 9.217822 9.031900 9.129272 8.795796
## ENSG00000000460 7.402051 8.392672 7.177641 8.066668 8.656107
## ENSG00000000938 10.139142 8.430335 8.337217 8.964144 8.971527
## ENSG00000000971 12.144935 10.191008 11.952959 9.673892 9.778397
## ENSG00000001036 10.341637 10.420020 10.289594 10.281137 10.391436
## ENSG00000001084 9.295030 8.900146 9.200154 9.180694 9.221930
## ENSG00000001167 8.899336 9.471624 9.227674 9.078129 10.028288
## SRR10912147 SRR10912148 SRR10912149 SRR10912150 SRR10912151
## ENSG00000000003 10.984098 11.419437 10.896370 9.860183 9.919823
## ENSG00000000005 6.700062 6.226983 5.657905 5.657905 9.117115
## ENSG00000000419 10.679558 10.030078 10.530196 9.784860 10.341940
## ENSG00000000457 8.633689 8.824891 8.475233 8.598698 8.898314
## ENSG00000000460 7.395527 8.772661 7.720103 7.824260 6.559439
## ENSG00000000938 6.700062 7.824077 8.083384 7.539390 8.799438
## ENSG00000000971 9.338240 8.928647 11.302693 11.403038 12.637810
## ENSG00000001036 9.715562 10.989807 10.136391 9.704060 10.482927
## ENSG00000001084 9.067209 9.228072 9.115244 8.555554 9.745390
## ENSG00000001167 8.152816 8.548046 7.949802 8.271179 8.668850
## SRR10912152 SRR10912153 SRR10912154 SRR10912155 SRR10912156
## ENSG00000000003 10.078535 9.852815 10.388402 10.079501 10.611783
## ENSG00000000005 5.657905 5.657905 5.657905 5.911852 5.657905
## ENSG00000000419 10.259810 9.957853 10.290172 10.147716 9.596196
## ENSG00000000457 8.836992 8.558326 9.535682 8.058282 8.434271
## ENSG00000000460 8.548649 8.248089 7.897485 7.859239 7.377141
## ENSG00000000938 7.454627 7.996913 8.775385 8.859813 9.795498
## ENSG00000000971 9.647801 8.380668 9.766278 10.540512 10.740995
## ENSG00000001036 10.034946 9.980854 10.135844 10.385985 9.278702
## ENSG00000001084 8.153185 8.961764 8.649668 9.824866 8.309312
## ENSG00000001167 8.299080 8.384380 8.473338 8.762909 8.366951
## SRR10912157 SRR10912158 SRR10912159 SRR10912160 SRR10912161
## ENSG00000000003 11.795582 11.389175 7.463260 8.449949 5.657905
## ENSG00000000005 5.907736 5.657905 5.657905 5.657905 5.657905
## ENSG00000000419 10.287970 10.563573 9.092711 9.122085 5.657905
## ENSG00000000457 8.777955 8.552365 6.973376 8.080911 5.657905
## ENSG00000000460 8.639941 8.412397 7.810745 7.763193 5.657905
## ENSG00000000938 6.290565 7.764174 11.197375 11.934724 11.921847
## ENSG00000000971 7.195523 8.474792 11.224753 10.904849 10.797183
## ENSG00000001036 10.550255 9.306829 10.587098 10.629538 10.797183
## ENSG00000001084 9.542411 9.288514 8.311946 8.905904 9.011239
## ENSG00000001167 8.875595 8.397589 8.830090 8.740883 5.657905
## SRR10912162 SRR10912163 SRR10912164 SRR10912165 SRR10912166
## ENSG00000000003 7.943794 10.470963 10.605714 10.772655 5.657905
## ENSG00000000005 5.657905 5.657905 5.657905 5.657905 5.657905
## ENSG00000000419 9.325488 10.749076 10.105470 9.926235 5.657905
## ENSG00000000457 8.664924 8.612039 7.967889 8.958865 5.657905
## ENSG00000000460 8.351213 7.970324 8.427629 8.506545 5.657905
## ENSG00000000938 10.214300 8.446406 7.254249 6.790146 11.211340
## ENSG00000000971 10.116392 8.700865 8.427629 8.212602 11.211340
## ENSG00000001036 10.305962 10.482355 9.651894 9.800036 11.211340
## ENSG00000001084 8.920954 9.674335 8.846012 8.745098 5.657905
## ENSG00000001167 8.664924 9.544415 9.331168 9.379211 5.657905
## SRR10912167 SRR10912168 SRR10912169 SRR10912170 SRR10912171
## ENSG00000000003 8.180514 10.547861 8.892073 10.261833 10.917607
## ENSG00000000005 5.657905 5.657905 5.657905 5.833176 6.128923
## ENSG00000000419 8.975025 9.892837 10.827876 9.939263 10.724170
## ENSG00000000457 8.521530 8.840929 8.410319 8.561306 8.410402
## ENSG00000000460 7.517007 7.740980 8.389061 8.253949 8.627841
## ENSG00000000938 12.081720 8.816368 8.525899 7.856368 6.977212
## ENSG00000000971 11.537685 10.113472 8.329650 9.293253 9.651555
## ENSG00000001036 10.560302 9.520768 9.799506 10.599937 10.217742
## ENSG00000001084 8.597257 8.543741 9.788173 9.166928 9.180999
## ENSG00000001167 9.195511 9.555541 8.323244 9.571613 8.751902
## SRR10912172
## ENSG00000000003 8.903937
## ENSG00000000005 5.657905
## ENSG00000000419 9.165333
## ENSG00000000457 8.142573
## ENSG00000000460 7.627294
## ENSG00000000938 11.488695
## ENSG00000000971 10.114272
## ENSG00000001036 10.389456
## ENSG00000001084 9.028052
## ENSG00000001167 8.914531meanSdPlot(assay(ntd))
meanSdPlot(assay(vsd))
#meanSdPlot(assay(rld))
vsd <- vst(dds_Treatment, blind=FALSE)
#rld <- rlog(dds_Tissue, blind=FALSE)
ntd <- normTransform(dds_Treatment)
library("vsn")
head(assay(vsd), 10)## SRR10912052 SRR10912053 SRR10912054 SRR10912055 SRR10912056
## ENSG00000000003 9.635131 11.394629 10.343799 10.840225 9.761264
## ENSG00000000005 5.198234 7.026823 7.118182 5.198234 5.198234
## ENSG00000000419 10.109271 9.921134 10.126099 10.455749 10.168354
## ENSG00000000457 8.693937 8.635247 8.612133 8.486984 8.880130
## ENSG00000000460 8.252147 7.544796 8.188231 7.701604 7.817548
## ENSG00000000938 8.695669 8.786432 8.099575 6.630792 8.589261
## ENSG00000000971 11.331935 12.173337 12.112666 9.624793 13.724290
## ENSG00000001036 10.294303 10.741870 10.512750 9.686202 10.808320
## ENSG00000001084 8.729866 9.764873 9.342162 9.374682 9.042639
## ENSG00000001167 9.746450 9.089312 9.137673 9.260998 8.992501
## SRR10912057 SRR10912058 SRR10912059 SRR10912060 SRR10912061
## ENSG00000000003 10.659093 9.044602 8.358228 10.615424 10.169648
## ENSG00000000005 5.405480 5.667645 5.198234 6.463288 5.198234
## ENSG00000000419 9.986798 9.748198 9.221749 10.339659 9.772742
## ENSG00000000457 8.912778 8.655996 8.247392 8.658319 8.616751
## ENSG00000000460 8.622284 7.414579 7.854441 6.963492 5.198234
## ENSG00000000938 7.477206 10.038424 10.816167 7.702903 8.572418
## ENSG00000000971 11.178017 13.157804 11.131692 12.206168 11.399017
## ENSG00000001036 10.417759 10.905977 11.136015 10.396621 9.990700
## ENSG00000001084 9.622296 9.676591 9.593633 8.274971 8.540538
## ENSG00000001167 9.802621 9.046753 8.500656 9.235261 8.928968
## SRR10912062 SRR10912063 SRR10912064 SRR10912065 SRR10912066
## ENSG00000000003 10.451730 9.650686 11.388246 10.611081 10.638724
## ENSG00000000005 5.198234 6.227182 5.198234 5.655517 5.650143
## ENSG00000000419 9.960810 9.541530 9.956935 9.618100 10.471560
## ENSG00000000457 7.358650 7.926039 9.132185 9.262732 8.257383
## ENSG00000000460 7.989102 7.907240 8.620390 8.126036 8.481714
## ENSG00000000938 8.775771 9.444829 8.125104 8.359504 6.611403
## ENSG00000000971 12.163049 12.326997 12.137347 12.303419 10.757039
## ENSG00000001036 10.912119 10.636659 10.732071 10.479376 10.401573
## ENSG00000001084 9.382642 10.341621 9.580560 9.533539 9.123281
## ENSG00000001167 8.994202 9.016373 9.183171 9.079016 9.486041
## SRR10912067 SRR10912068 SRR10912069 SRR10912070 SRR10912071
## ENSG00000000003 9.462146 9.703693 8.852233 10.630851 9.117855
## ENSG00000000005 5.672231 5.627259 5.637854 5.489889 5.559554
## ENSG00000000419 10.186043 10.055273 9.684774 10.313201 10.287538
## ENSG00000000457 8.775580 8.799971 8.880727 8.728089 8.423869
## ENSG00000000460 7.706897 7.888726 7.631715 8.398392 7.633631
## ENSG00000000938 11.248619 8.927340 10.104017 7.482233 9.408296
## ENSG00000000971 10.323069 11.305849 11.552952 11.363773 13.649729
## ENSG00000001036 10.466660 10.326613 10.295137 10.273357 10.240420
## ENSG00000001084 10.041459 8.970543 9.366316 9.389353 9.095545
## ENSG00000001167 9.046416 9.504868 9.229408 9.082612 8.860552
## SRR10912072 SRR10912073 SRR10912074 SRR10912075 SRR10912076
## ENSG00000000003 9.746388 11.794152 10.340581 9.477069 9.311941
## ENSG00000000005 5.434435 6.364305 7.146107 5.198234 6.998363
## ENSG00000000419 9.828634 10.377110 10.326082 9.819279 9.919435
## ENSG00000000457 8.502399 9.194608 9.069602 9.502416 8.963036
## ENSG00000000460 8.399638 8.417694 8.786101 8.748361 7.765716
## ENSG00000000938 8.601892 6.619047 8.848657 8.635167 9.086898
## ENSG00000000971 12.701249 8.158492 10.966005 10.336551 12.355547
## ENSG00000001036 10.287702 9.697241 10.340581 9.947715 9.917639
## ENSG00000001084 9.085828 9.621819 9.734599 9.591704 9.702874
## ENSG00000001167 9.142562 10.406904 9.837579 9.583812 8.970013
## SRR10912077 SRR10912078 SRR10912079 SRR10912080 SRR10912081
## ENSG00000000003 9.174712 11.674195 11.640170 10.929598 9.844063
## ENSG00000000005 5.989768 7.633864 6.196801 5.198234 5.427470
## ENSG00000000419 9.649499 10.478871 10.369345 11.678990 10.821512
## ENSG00000000457 8.778541 8.606720 8.638979 8.584815 8.632974
## ENSG00000000460 7.758695 8.360635 7.755565 8.394797 7.633267
## ENSG00000000938 9.684306 7.526275 6.997363 7.681926 9.587321
## ENSG00000000971 12.550873 11.240292 12.566202 11.298481 12.403271
## ENSG00000001036 10.474689 10.471435 10.259963 11.839145 11.372574
## ENSG00000001084 9.365410 8.927061 9.276856 10.331155 9.830859
## ENSG00000001167 9.053531 10.342846 9.637218 9.466062 8.796789
## SRR10912082 SRR10912083 SRR10912084 SRR10912085 SRR10912086
## ENSG00000000003 10.708167 10.601412 10.362631 9.814877 10.027085
## ENSG00000000005 5.198234 5.198234 5.198234 8.198750 9.849574
## ENSG00000000419 10.071091 10.431703 10.286364 10.246700 9.892890
## ENSG00000000457 9.417572 8.627173 8.716299 8.884176 8.608305
## ENSG00000000460 9.262932 8.061438 8.595586 7.740423 7.695662
## ENSG00000000938 8.110680 7.762126 8.653011 8.913911 8.990321
## ENSG00000000971 11.361986 12.602566 11.588347 13.275553 12.934971
## ENSG00000001036 9.464097 9.982106 10.527061 10.428176 10.602164
## ENSG00000001084 9.942468 9.319070 9.801830 9.694404 9.704009
## ENSG00000001167 10.097851 9.941472 9.930761 9.786056 8.786664
## SRR10912087 SRR10912088 SRR10912089 SRR10912090 SRR10912091
## ENSG00000000003 10.086073 10.607756 10.728015 10.088830 9.812958
## ENSG00000000005 8.206353 6.289888 5.661170 5.410677 6.277873
## ENSG00000000419 10.136875 10.026697 9.730718 9.707086 9.513685
## ENSG00000000457 9.044031 9.334401 8.923268 8.539413 8.627809
## ENSG00000000460 7.366673 8.439982 9.026557 7.464770 7.528132
## ENSG00000000938 10.231003 7.401353 7.310425 7.712579 7.548205
## ENSG00000000971 12.340936 10.098152 9.148412 11.320341 12.666925
## ENSG00000001036 10.304620 10.089411 11.188008 9.982390 10.057176
## ENSG00000001084 10.016813 9.508697 9.417403 8.908539 8.571865
## ENSG00000001167 9.182355 9.331927 9.376486 8.803619 8.687485
## SRR10912092 SRR10912093 SRR10912094 SRR10912095 SRR10912096
## ENSG00000000003 9.644322 9.729994 9.529662 10.067375 10.563996
## ENSG00000000005 8.501919 5.198234 5.468650 5.492601 5.560411
## ENSG00000000419 9.402949 10.020641 9.616824 10.167077 9.583315
## ENSG00000000457 8.409479 9.390107 8.725896 8.659840 8.743066
## ENSG00000000460 7.488081 8.636475 8.231231 7.851354 8.269737
## ENSG00000000938 9.157826 6.886545 7.823975 7.358813 7.620814
## ENSG00000000971 12.883026 9.959229 8.513068 11.329637 10.781780
## ENSG00000001036 10.513337 9.942560 9.908991 10.204163 9.923921
## ENSG00000001084 9.513718 8.441743 9.021603 9.178779 9.461364
## ENSG00000001167 9.198821 9.081105 8.607872 9.481715 9.547805
## SRR10912097 SRR10912098 SRR10912099 SRR10912100 SRR10912101
## ENSG00000000003 11.907497 10.704689 11.757826 10.707248 10.093059
## ENSG00000000005 6.497883 5.412206 5.198234 5.198234 5.198234
## ENSG00000000419 9.415738 9.296297 10.146167 9.936246 9.487264
## ENSG00000000457 8.810950 8.550167 9.117088 8.805634 8.112872
## ENSG00000000460 7.268165 7.512139 8.886472 8.802646 7.286579
## ENSG00000000938 6.728786 7.692652 6.138172 7.395948 7.319552
## ENSG00000000971 12.429678 10.586625 7.356246 10.575134 11.518326
## ENSG00000001036 11.028457 9.285077 10.726300 9.560203 9.759156
## ENSG00000001084 7.880227 8.614699 10.015858 9.088961 8.966030
## ENSG00000001167 9.402698 8.955011 9.991885 9.327381 9.752020
## SRR10912102 SRR10912103 SRR10912104 SRR10912105 SRR10912106
## ENSG00000000003 9.331289 10.746147 10.638844 10.413793 10.017774
## ENSG00000000005 7.876887 5.479091 5.978438 5.498280 5.930750
## ENSG00000000419 8.977935 11.017936 10.478901 9.577689 9.536164
## ENSG00000000457 8.521122 8.816297 8.646889 8.825747 8.510412
## ENSG00000000460 7.025808 9.188256 8.361786 7.841170 8.233742
## ENSG00000000938 9.293018 7.104418 8.177590 7.694662 7.650235
## ENSG00000000971 12.575053 8.897150 10.085501 9.060372 10.627789
## ENSG00000001036 10.184019 10.858370 9.986635 10.232489 9.819448
## ENSG00000001084 9.825562 9.580005 9.629273 9.993754 9.976774
## ENSG00000001167 8.905632 9.593131 9.857808 9.825046 9.410763
## SRR10912107 SRR10912108 SRR10912109 SRR10912110 SRR10912111
## ENSG00000000003 9.796485 9.961830 11.041643 10.937511 9.296541
## ENSG00000000005 6.062247 5.557402 6.032298 6.414506 5.572553
## ENSG00000000419 9.639392 9.505545 9.363771 9.541241 10.347239
## ENSG00000000457 8.059507 9.088547 8.532219 8.821393 8.891409
## ENSG00000000460 7.797877 8.075244 7.533074 7.908898 8.643070
## ENSG00000000938 8.469702 7.811697 7.384392 6.768681 7.600389
## ENSG00000000971 8.863095 10.041583 10.467734 9.738373 9.242120
## ENSG00000001036 9.616149 10.266742 10.593880 9.679420 10.260181
## ENSG00000001084 9.185838 9.263157 9.106418 10.005443 8.744861
## ENSG00000001167 8.884087 9.217988 9.321587 9.477158 9.919501
## SRR10912112 SRR10912113 SRR10912114 SRR10912115 SRR10912116
## ENSG00000000003 10.523119 10.135530 9.958470 11.047404 10.881853
## ENSG00000000005 5.198234 5.980160 5.640500 5.884616 5.198234
## ENSG00000000419 9.598331 10.127527 10.336406 9.982054 9.899575
## ENSG00000000457 8.983249 8.767242 8.501820 8.882942 8.796903
## ENSG00000000460 8.820641 8.395432 8.107402 8.355736 7.721291
## ENSG00000000938 7.240943 7.831319 7.162926 7.219798 6.867067
## ENSG00000000971 9.542843 10.274887 9.889170 10.136728 10.761600
## ENSG00000001036 10.020973 10.416606 10.089639 9.431178 9.782956
## ENSG00000001084 9.038661 9.265641 9.484285 8.455946 8.858263
## ENSG00000001167 9.167469 10.002779 10.110687 9.958681 9.294703
## SRR10912117 SRR10912118 SRR10912119 SRR10912120 SRR10912121
## ENSG00000000003 10.849212 9.747802 10.131316 10.544054 8.942871
## ENSG00000000005 5.432279 8.198908 5.716403 5.198234 5.198234
## ENSG00000000419 10.114173 8.899338 9.744689 10.006138 11.117377
## ENSG00000000457 8.875014 8.727325 8.514078 8.614346 8.469639
## ENSG00000000460 8.538863 7.540210 8.411184 8.177136 8.244182
## ENSG00000000938 7.225280 8.922969 8.053026 7.567644 7.953960
## ENSG00000000971 9.720404 11.867829 11.313454 11.952069 10.322610
## ENSG00000001036 10.396312 9.924899 10.061858 10.243250 10.165468
## ENSG00000001084 9.927132 9.497965 9.832294 8.922145 9.360638
## ENSG00000001167 9.067933 8.820972 10.446714 9.585889 8.799974
## SRR10912122 SRR10912123 SRR10912124 SRR10912125 SRR10912126
## ENSG00000000003 9.995809 10.051583 9.520825 10.747392 9.960777
## ENSG00000000005 5.502618 5.419802 5.198234 5.966735 5.696518
## ENSG00000000419 9.749623 10.295417 9.520825 9.659084 9.676969
## ENSG00000000457 8.322764 8.493248 8.131113 8.724467 8.504320
## ENSG00000000460 8.121784 8.561764 7.065470 8.238870 7.644770
## ENSG00000000938 7.511783 6.725268 10.670943 7.506350 7.940204
## ENSG00000000971 11.144323 11.379225 11.992983 11.339232 11.382685
## ENSG00000001036 10.515201 10.321106 10.259540 10.697937 10.106021
## ENSG00000001084 9.227279 9.267072 9.714842 9.356949 8.595980
## ENSG00000001167 9.537846 9.357919 9.781923 9.460719 8.642364
## SRR10912127 SRR10912128 SRR10912129 SRR10912130 SRR10912131
## ENSG00000000003 11.534463 11.409715 10.489673 11.186739 11.163333
## ENSG00000000005 6.800841 5.976597 7.423524 6.213450 6.240671
## ENSG00000000419 9.655482 9.490512 9.466251 9.584142 9.790377
## ENSG00000000457 9.046576 8.813700 8.595021 8.865231 8.730175
## ENSG00000000460 7.118566 7.615514 6.930873 7.157304 7.052479
## ENSG00000000938 6.492453 6.961793 7.450526 6.736734 6.297623
## ENSG00000000971 12.259767 12.690303 13.114110 12.076900 13.142028
## ENSG00000001036 10.669254 10.438879 10.682945 10.604614 10.585067
## ENSG00000001084 9.036354 9.640207 8.336611 8.519058 8.634317
## ENSG00000001167 9.580286 9.540984 9.240283 9.230211 9.616991
## SRR10912132 SRR10912133 SRR10912134 SRR10912135 SRR10912136
## ENSG00000000003 11.070523 11.337754 11.548295 11.145485 11.062804
## ENSG00000000005 8.119712 5.994681 6.238808 6.386620 7.021958
## ENSG00000000419 9.577299 9.759631 9.785589 9.725496 9.653751
## ENSG00000000457 8.730744 8.943014 9.069817 9.036626 8.700004
## ENSG00000000460 7.267416 7.321361 7.502308 7.290197 6.994754
## ENSG00000000938 7.879936 6.601457 6.071668 6.626184 6.269433
## ENSG00000000971 12.672010 13.094924 12.711257 12.570939 12.893964
## ENSG00000001036 10.319821 10.065583 10.559589 10.654218 10.536619
## ENSG00000001084 8.900332 8.902465 9.222710 8.616691 8.434975
## ENSG00000001167 9.371456 9.595831 9.758357 9.696409 9.421845
## SRR10912137 SRR10912138 SRR10912139 SRR10912140 SRR10912141
## ENSG00000000003 11.187095 10.435172 10.532387 9.294067 10.051733
## ENSG00000000005 7.024070 5.198234 5.454690 5.380794 5.514578
## ENSG00000000419 9.567217 10.600018 10.214951 9.704969 10.322469
## ENSG00000000457 8.786111 8.182587 8.482983 8.068520 8.309696
## ENSG00000000460 7.168708 7.743096 8.145190 7.506874 7.916058
## ENSG00000000938 6.270432 9.027989 8.892413 9.811852 9.459271
## ENSG00000000971 12.709759 11.458503 9.709828 10.993011 9.999338
## ENSG00000001036 10.448255 10.867579 10.398570 10.844987 9.612834
## ENSG00000001084 8.207451 8.697029 9.248947 9.517370 9.396611
## ENSG00000001167 9.739384 9.381255 9.260761 9.113199 8.937856
## SRR10912142 SRR10912143 SRR10912144 SRR10912145 SRR10912146
## ENSG00000000003 9.433316 11.408671 10.793277 9.873178 10.403139
## ENSG00000000005 6.704740 5.638440 5.198234 5.279740 5.698854
## ENSG00000000419 9.347726 10.051484 9.711179 9.516300 10.013252
## ENSG00000000457 8.525830 9.154332 8.960173 9.061977 8.712166
## ENSG00000000460 7.204769 8.284444 6.953668 7.933998 8.564600
## ENSG00000000938 10.104815 8.324655 8.225137 8.889174 8.896917
## ENSG00000000971 12.136216 10.157862 11.943009 9.626976 9.734648
## ENSG00000001036 10.311701 10.391633 10.258584 10.249950 10.362494
## ENSG00000001084 9.234691 8.821987 9.135921 9.115633 9.158611
## ENSG00000001167 8.821136 9.417943 9.164594 9.008539 9.991297
## SRR10912147 SRR10912148 SRR10912149 SRR10912150 SRR10912151
## ENSG00000000003 10.964761 11.405082 10.875840 9.818768 9.880033
## ENSG00000000005 6.411086 5.864005 5.198234 5.198234 9.049280
## ENSG00000000419 10.655759 9.993132 10.503854 9.741300 10.312011
## ENSG00000000457 8.540857 8.742822 8.372521 8.503761 8.820062
## ENSG00000000460 7.197505 8.687772 7.556365 7.670443 6.249347
## ENSG00000000938 6.411086 7.670243 7.952076 7.357203 8.716005
## ENSG00000000971 9.279599 8.851924 11.287143 11.388521 12.631602
## ENSG00000001036 9.669934 10.970546 10.102000 9.658081 10.455726
## ENSG00000001084 8.997119 9.165009 9.047326 8.457963 9.700663
## ENSG00000001167 8.027035 8.449987 7.807268 8.154352 8.578090
## SRR10912152 SRR10912153 SRR10912154 SRR10912155 SRR10912156
## ENSG00000000003 10.042776 9.811195 10.359400 10.043765 10.586862
## ENSG00000000005 5.198234 5.198234 5.198234 5.495895 5.198234
## ENSG00000000419 10.228170 9.919067 10.259175 10.113586 9.546784
## ENSG00000000457 8.755564 8.460908 9.484239 7.924925 8.328855
## ENSG00000000460 8.450627 8.129560 7.750337 7.708638 7.177021
## ENSG00000000938 7.263234 7.858429 8.690645 8.779580 9.752248
## ENSG00000000971 9.600061 8.271616 9.722173 10.514354 10.718171
## ENSG00000001036 9.998121 9.942663 10.101440 10.356935 9.217709
## ENSG00000001084 8.027433 8.886677 8.557783 9.782458 8.195246
## ENSG00000001167 8.184279 8.275583 8.370502 8.677484 8.256950
## SRR10912157 SRR10912158 SRR10912159 SRR10912160 SRR10912161
## ENSG00000000003 11.784494 11.374519 7.272821 8.345575 5.198234
## ENSG00000000005 5.491074 5.198234 5.198234 5.198234 5.198234
## ENSG00000000419 10.256927 10.537821 9.023783 9.054471 5.198234
## ENSG00000000457 8.693355 8.454576 6.722934 7.949402 5.198234
## ENSG00000000460 8.547480 8.305510 7.655669 7.603622 5.198234
## ENSG00000000938 5.937985 7.604696 11.180662 11.924648 11.911681
## ENSG00000000971 6.973770 8.372051 11.208350 10.884438 10.775216
## ENSG00000001036 10.524269 9.246958 10.561754 10.604915 10.775216
## ENSG00000001084 9.491197 9.227914 8.198070 8.828037 8.938538
## ENSG00000001167 8.796179 8.289696 8.748296 8.654235 5.198234
## SRR10912162 SRR10912163 SRR10912164 SRR10912165 SRR10912166
## ENSG00000000003 7.800737 10.443540 10.580689 10.750318 5.198234
## ENSG00000000005 5.198234 5.198234 5.198234 5.198234 5.198234
## ENSG00000000419 9.266350 10.726377 10.070354 9.886616 5.198234
## ENSG00000000457 8.573934 8.517910 7.826923 8.883636 5.198234
## ENSG00000000460 8.240116 7.829567 8.321768 8.405857 5.198234
## ENSG00000000938 10.181671 8.341797 7.039671 6.514250 11.194786
## ENSG00000000971 10.081534 8.611954 8.321768 8.091417 11.194786
## ENSG00000001036 10.275295 10.455143 9.604284 9.756917 11.194786
## ENSG00000001084 8.843845 9.627433 8.765058 8.658684 5.198234
## ENSG00000001167 8.573934 9.493270 9.272253 9.322139 5.198234
## SRR10912167 SRR10912168 SRR10912169 SRR10912170 SRR10912171
## ENSG00000000003 8.056880 10.521832 8.813503 10.230236 10.897372
## ENSG00000000005 5.198234 5.198234 5.198234 5.403728 5.749691
## ENSG00000000419 8.900585 9.852319 10.806364 9.899989 10.701083
## ENSG00000000457 8.421799 8.759707 8.303292 8.464073 8.303381
## ENSG00000000460 7.332423 7.579272 8.280586 8.135854 8.534660
## ENSG00000000938 12.072613 8.733844 8.426445 7.705505 6.727285
## ENSG00000000971 11.524448 10.078545 8.217034 9.232842 9.603934
## ENSG00000001036 10.534493 9.468812 9.756371 10.574814 10.185189
## ENSG00000001084 8.502233 8.445411 9.744710 9.101276 9.115951
## ENSG00000001167 9.131081 9.504773 8.210174 9.521386 8.665868
## SRR10912172
## ENSG00000000003 8.825971
## ENSG00000000005 5.198234
## ENSG00000000419 9.099612
## ENSG00000000457 8.015989
## ENSG00000000460 7.454280
## ENSG00000000938 11.475006
## ENSG00000000971 10.079364
## ENSG00000001036 10.360475
## ENSG00000001084 8.956144
## ENSG00000001167 8.837100meanSdPlot(assay(ntd))
meanSdPlot(assay(vsd))
#meanSdPlot(assay(rld))Heatmaps
library("pheatmap")
select <- order(rowMeans(counts(dds_Tissue,normalized=TRUE)),
decreasing=TRUE)[1:20]
df <- as.data.frame(colData(dds_Tissue)[,c("Tissue","Treatment")])
pheatmap(assay(ntd)[select,], cluster_rows=TRUE, show_rownames=FALSE,
cluster_cols=TRUE, annotation_col=df)
pheatmap(assay(vsd)[select,], cluster_rows=TRUE, show_rownames=FALSE,
cluster_cols=TRUE, annotation_col=df)
#pheatmap(assay(rld)[select,], cluster_rows=FALSE, show_rownames=FALSE,
# cluster_cols=FALSE, annotation_col=df)sampleDists <- dist(t(assay(vsd)))
library("RColorBrewer")
sampleDistMatrix <- as.matrix(sampleDists)
rownames(sampleDistMatrix) <- paste(vsd$condition, vsd$type, sep="-")
colnames(sampleDistMatrix) <- NULL
colors <- colorRampPalette( rev(brewer.pal(9, "Blues")) )(255)
pheatmap(sampleDistMatrix,
clustering_distance_rows=sampleDists,
clustering_distance_cols=sampleDists,
col=colors)
PCA
plotPCA(vsd, intgroup=c("Tissue"))
if (!require("BiocManager", quietly = TRUE))
install.packages("BiocManager")
BiocManager::install("clusterProfiler")## Bioconductor version 3.17 (BiocManager 1.30.21), R 4.3.0 (2023-04-21)## Warning: package(s) not installed when version(s) same as or greater than current; use
## `force = TRUE` to re-install: 'clusterProfiler'## Old packages: 'DESeq2', 'digest', 'downlit', 'GenomeInfoDb', 'GenomicFeatures',
## 'IRanges', 'jsonlite', 'MatrixGenerics', 'pkgbuild'Random
results(dds_Tissue, contrast=c("Tissue","Serous_EOC","Benign_Tissue"))## log2 fold change (MLE): Tissue Serous_EOC vs Benign_Tissue
## Wald test p-value: Tissue Serous EOC vs Benign Tissue
## DataFrame with 19975 rows and 6 columns
## baseMean log2FoldChange lfcSE stat pvalue
## <numeric> <numeric> <numeric> <numeric> <numeric>
## ENSG00000000003 1359.4263 -0.790170 0.248072 -3.185252 1.44628e-03
## ENSG00000000005 28.6264 -0.512972 0.913941 -0.561275 5.74610e-01
## ENSG00000000419 945.1298 0.408815 0.156660 2.609561 9.06584e-03
## ENSG00000000457 339.0323 -0.139602 0.124784 -1.118752 2.63246e-01
## ENSG00000000460 189.3951 1.346737 0.231249 5.823762 5.75375e-09
## ... ... ... ... ... ...
## ENSG00000273294 18.6100019 0.115738 0.774778 0.1493822 0.881252
## ENSG00000273331 1.6135503 2.457218 1.142924 2.1499398 0.031560
## ENSG00000273398 5.7510137 0.562653 0.536606 1.0485403 0.294390
## ENSG00000273439 22.4666648 0.669630 0.415259 1.6125583 0.106841
## ENSG00000273452 0.0657423 0.197680 4.954781 0.0398967 0.968175
## padj
## <numeric>
## ENSG00000000003 3.91623e-03
## ENSG00000000005 6.64751e-01
## ENSG00000000419 1.99003e-02
## ENSG00000000457 3.58693e-01
## ENSG00000000460 5.14674e-08
## ... ...
## ENSG00000273294 0.9141591
## ENSG00000273331 0.0593248
## ENSG00000273398 0.3925197
## ENSG00000273439 0.1682453
## ENSG00000273452 NAMore LFC Shrinkage
resApeT <- lfcShrink(dds_Tissue, coef=2, type="apeglm", lfcThreshold=1)## using 'apeglm' for LFC shrinkage. If used in published research, please cite:
## Zhu, A., Ibrahim, J.G., Love, M.I. (2018) Heavy-tailed prior distributions for
## sequence count data: removing the noise and preserving large differences.
## Bioinformatics. https://doi.org/10.1093/bioinformatics/bty895## computing FSOS 'false sign or small' s-values (T=1)plotMA(resApeT, ylim=c(-3,3), cex=.8)## thresholding s-values on alpha=0.005 to color pointsabline(h=c(-1,1), col="dodgerblue", lwd=2)
resAshT <- lfcShrink(dds_Tissue, coef=2, type="ashr", lfcThreshold=1)## using 'ashr' for LFC shrinkage. If used in published research, please cite:
## Stephens, M. (2016) False discovery rates: a new deal. Biostatistics, 18:2.
## https://doi.org/10.1093/biostatistics/kxw041## computing FSOS 'false sign or small' s-values (T=1)plotMA(resAshT, ylim=c(-3,3), cex=.8)## thresholding s-values on alpha=0.005 to color pointsabline(h=c(-1,1), col="dodgerblue", lwd=2)
Boxplot
par(mar=c(8,5,2,2))
boxplot(log10(assays(dds_Tissue)[["cooks"]]), range=0, las=2)
Dispersion Plot
plotDispEsts(dds_Tissue)
Other
metadata(resLFC_Serous_EOC)$alpha## [1] 0.1metadata(resLFC_Serous_EOC)$filterThreshold## 3.896759%
## 0.2997622plot(metadata(resLFC_Serous_EOC)$filterNumRej,
type="b", ylab="number of rejections",
xlab="quantiles of filter")
lines(metadata(resLFC_Serous_EOC)$lo.fit, col="red")
abline(v=metadata(resLFC_Serous_EOC)$filterTheta)
resNoFilt <- results(dds_Tissue, independentFiltering=FALSE)
addmargins(table(filtering=(resLFC_Serous_EOC$padj < .1),
noFiltering=(resNoFilt$padj < .1)))## noFiltering
## filtering FALSE TRUE Sum
## FALSE 8060 0 8060
## TRUE 84 11052 11136
## Sum 8144 11052 19196metadata(resLFC_Ascites)$alpha## [1] 0.1metadata(resLFC_Ascites)$filterThreshold## 0.02002503%
## 0.005648919plot(metadata(resLFC_Ascites)$filterNumRej,
type="b", ylab="number of rejections",
xlab="quantiles of filter")
lines(metadata(resLFC_Ascites)$lo.fit, col="red")
abline(v=metadata(resLFC_Ascites)$filterTheta)
resNoFilt <- results(dds_Tissue, independentFiltering=FALSE)
addmargins(table(filtering=(resLFC_Ascites$padj < .1),
noFiltering=(resNoFilt$padj < .1)))## noFiltering
## filtering FALSE TRUE Sum
## FALSE 5118 2094 7212
## TRUE 3801 8958 12759
## Sum 8919 11052 19971par(mfrow=c(2,2),mar=c(2,2,1,1))
ylim <- c(-2.5,2.5)
resGA <- results(dds_Tissue, lfcThreshold=.5, altHypothesis="greaterAbs")
resLA <- results(dds_Tissue, lfcThreshold=.5, altHypothesis="lessAbs")
resG <- results(dds_Tissue, lfcThreshold=.5, altHypothesis="greater")
resL <- results(dds_Tissue, lfcThreshold=.5, altHypothesis="less")
drawLines <- function() abline(h=c(-.5,.5),col="dodgerblue",lwd=2)
plotMA(resGA, ylim=ylim); drawLines()
plotMA(resLA, ylim=ylim); drawLines()
plotMA(resG, ylim=ylim); drawLines()
plotMA(resL, ylim=ylim); drawLines()
#mcols(dds_Tissue,use.names=TRUE)[1:4,1:4]
#substr(names(mcols(dds_Tissue)),1,10)
#mcols(mcols(dds_Tissue), use.names=TRUE)[1:4,]
#head(assays(dds_Tissue)[["mu"]])
#head(assays(dds_Tissue)[["cooks"]])
#head(dispersions(dds_Tissue))
#head(mcols(dds_Tissue)$dispersion)
#sizeFactors(dds_Tissue)
#head(coef(dds_Tissue))
#attr(dds_Tissue, "betaPriorVar")
#priorInfo(resLFC)library("AnnotationDbi")
library("org.Hs.eg.db")## final_data_Serous <- read.csv('Final_Data_LFC_Serous_EOC.csv')
final_data_Serous$symbol <- mapIds(org.Hs.eg.db, keys = final_data_Serous$X, column = 'SYMBOL', keytype = 'ENSEMBL')## 'select()' returned 1:many mapping between keys and columnsfinal_data_Serous <- final_data_Serous[order(final_data_Serous$pvalue),]
head(final_data_Serous, 50)## X baseMean log2FoldChange lfcSE pvalue
## 2280 ENSG00000100380 5724.96499 -2.366535 0.1601225 6.747169e-50
## 9682 ENSG00000154545 61.40843 11.232181 3.5294452 9.463669e-50
## 2216 ENSG00000100227 1539.32309 -1.927756 0.1436224 2.460812e-41
## 7803 ENSG00000139734 173.22782 4.271019 0.3212016 9.455479e-41
## 1090 ENSG00000069966 566.39767 -2.153479 0.1634850 4.286271e-40
## 16378 ENSG00000196839 229.11553 2.378881 0.1868901 8.342579e-38
## 1770 ENSG00000088325 1169.30542 4.383619 0.3511531 1.148710e-37
## 7554 ENSG00000138180 416.12170 4.765852 0.3811779 2.852664e-37
## 18107 ENSG00000215298 33.08025 10.777277 4.0612932 6.886416e-37
## 12487 ENSG00000169679 391.34577 4.145891 0.3398087 6.577500e-36
## 9046 ENSG00000148773 1080.00392 4.341723 0.3616491 3.673715e-35
## 11477 ENSG00000165480 121.06319 3.984229 0.3236120 7.013135e-35
## 5110 ENSG00000121152 211.51483 4.233840 0.3543575 6.315500e-34
## 6423 ENSG00000131747 1159.48202 4.176065 0.3611308 8.299477e-33
## 12935 ENSG00000171700 461.50253 2.475897 0.2162880 2.153309e-31
## 5392 ENSG00000123843 129.27485 -3.469274 0.2995853 2.190238e-31
## 3687 ENSG00000109805 334.73959 4.094134 0.3615999 2.353862e-31
## 11271 ENSG00000164611 496.61379 3.671111 0.3229150 4.220011e-31
## 4009 ENSG00000112242 575.60515 2.399186 0.2111830 4.720339e-31
## 1844 ENSG00000089685 598.18916 4.314461 0.3844641 5.260845e-31
## 10002 ENSG00000157456 350.27085 3.846288 0.3409431 7.250601e-31
## 8275 ENSG00000143228 232.28054 3.839719 0.3421960 1.470819e-30
## 4081 ENSG00000112742 275.04632 4.183609 0.3760899 2.478286e-30
## 1372 ENSG00000077152 266.52244 3.316052 0.2964858 3.620661e-30
## 1197 ENSG00000072571 203.57910 3.767188 0.3388903 3.821498e-30
## 2420 ENSG00000101057 791.63169 4.677631 0.4245496 5.013545e-30
## 2516 ENSG00000101440 25.18471 -4.738041 0.4179163 5.225295e-30
## 4351 ENSG00000115163 159.78382 4.494506 0.4012504 9.625444e-30
## 923 ENSG00000065328 144.22937 3.772094 0.3401611 1.043422e-29
## 2523 ENSG00000101447 334.37476 4.109514 0.3763048 1.208282e-29
## 441 ENSG00000024526 225.58850 4.359077 0.3985200 2.076110e-29
## 14873 ENSG00000183287 422.64699 -3.630607 0.3253830 2.888063e-29
## 16917 ENSG00000198759 1127.80522 4.974515 0.4605761 3.667403e-29
## 14703 ENSG00000182481 2184.86118 2.733599 0.2508030 4.656406e-29
## 14757 ENSG00000182749 560.73077 -2.090683 0.1887148 6.188864e-29
## 3837 ENSG00000111206 780.85397 3.453181 0.3189707 6.999166e-29
## 16001 ENSG00000188730 126.42062 -5.363570 0.4823545 9.363982e-29
## 5377 ENSG00000123610 152.47400 5.081451 0.4642499 1.048219e-28
## 1198 ENSG00000072609 646.72286 1.419510 0.1293494 1.514072e-28
## 3846 ENSG00000111247 244.58649 3.034891 0.2809467 1.542488e-28
## 12834 ENSG00000171241 192.96019 3.138808 0.2910370 1.771562e-28
## 8938 ENSG00000147852 1134.30223 -3.043523 0.2772773 1.792106e-28
## 15017 ENSG00000183856 519.38832 4.345821 0.4069458 2.576997e-28
## 8223 ENSG00000142945 495.98642 3.918577 0.3673788 2.604511e-28
## 3141 ENSG00000105664 1740.27329 7.259892 0.6897367 2.897477e-28
## 18234 ENSG00000221955 679.22724 3.435748 0.3214557 3.241857e-28
## 12974 ENSG00000171848 565.55775 3.045456 0.2841631 3.339219e-28
## 9944 ENSG00000156970 295.96475 3.459288 0.3248491 4.514570e-28
## 5805 ENSG00000126787 267.17163 4.092292 0.3856053 9.048349e-28
## 8107 ENSG00000141905 1564.72016 -1.795886 0.1662445 1.107383e-27
## padj symbol
## 2280 9.083229e-46 ST13
## 9682 9.083229e-46 MAGED4
## 2216 1.574592e-37 POLDIP3
## 7803 4.537684e-37 DIAPH3
## 1090 1.645585e-36 GNB5
## 16378 2.669069e-34 ADA
## 1770 3.150091e-34 TPX2
## 7554 6.844966e-34 CEP55
## 18107 1.468796e-33 <NA>
## 12487 1.262617e-32 BUB1
## 9046 6.410967e-32 MKI67
## 11477 1.121868e-31 SKA3
## 5110 9.325564e-31 NCAPH
## 6423 1.137977e-29 TOP2A
## 12935 2.627739e-28 RGS19
## 5392 2.627739e-28 C4BPB
## 3687 2.657926e-28 NCAPG
## 11271 4.500407e-28 PTTG1
## 4009 4.769033e-28 E2F3
## 1844 5.049359e-28 BIRC5
## 10002 6.627740e-28 CCNB2
## 8275 1.283356e-27 NUF2
## 4081 2.068399e-27 TTK
## 1372 2.895925e-27 UBE2T
## 1197 2.934299e-27 HMMR
## 2420 3.701539e-27 MYBL2
## 2516 3.714991e-27 ASIP
## 4351 6.598929e-27 CENPA
## 923 6.906735e-27 MCM10
## 2523 7.731394e-27 FAM83D
## 441 1.285581e-26 DEPDC1
## 14873 1.732477e-26 CCBE1
## 16917 2.133317e-26 EGFL6
## 14703 2.628952e-26 KPNA2
## 14757 3.394327e-26 PAQR7
## 3837 3.732111e-26 FOXM1
## 16001 4.858135e-26 VWC2
## 5377 5.295163e-26 TNFAIP6
## 1198 7.402399e-26 CHFR
## 3846 7.402399e-26 RAD51AP1
## 12834 8.190776e-26 SHCBP1
## 8938 8.190776e-26 VLDLR
## 15017 1.136277e-25 IQGAP3
## 8223 1.136277e-25 KIF2C
## 3141 1.235999e-25 COMP
## 18234 1.352841e-25 SLC12A8
## 12974 1.363822e-25 RRM2
## 9944 1.805452e-25 BUB1B
## 5805 3.544737e-25 DLGAP5
## 8107 4.251464e-25 NFICfinal_data_Ascites <- read.csv('Final_Data_LFC_Ascites.csv')
final_data_Ascites$symbol <- mapIds(org.Hs.eg.db, keys = final_data_Ascites$X, column = 'SYMBOL', keytype = 'ENSEMBL')## 'select()' returned 1:many mapping between keys and columnsfinal_data_Ascites <- final_data_Ascites[order(final_data_Ascites$pvalue),]
head(final_data_Ascites, 50)## X baseMean log2FoldChange lfcSE pvalue
## 4410 ENSG00000115461 16817.08019 -6.766741 0.4209936 4.156538e-58
## 2280 ENSG00000100380 5724.96499 -2.749551 0.1781409 2.862340e-54
## 1209 ENSG00000072840 769.00105 -3.363033 0.2232460 7.753422e-52
## 2216 ENSG00000100227 1539.32309 -2.416056 0.1613224 2.936158e-51
## 4194 ENSG00000113658 1041.52626 -2.286800 0.1585170 1.053932e-47
## 4809 ENSG00000118640 2766.78334 3.871327 0.2781654 2.250980e-45
## 11261 ENSG00000164574 1551.57933 -3.680226 0.2639437 5.407182e-45
## 1924 ENSG00000091436 1251.49413 -3.798562 0.2736679 2.292355e-44
## 14389 ENSG00000180357 819.64498 -2.095509 0.1508520 3.046458e-44
## 1090 ENSG00000069966 566.39767 -2.530313 0.1842996 2.728464e-43
## 8092 ENSG00000141720 1522.08596 -2.120037 0.1546516 3.372399e-43
## 3841 ENSG00000111229 3753.62933 1.663285 0.1206307 4.872496e-43
## 7918 ENSG00000140543 138.51741 -2.592467 0.1912299 1.532077e-42
## 11027 ENSG00000163820 798.28237 -2.786760 0.2085839 1.632047e-41
## 15073 ENSG00000184207 578.75359 2.887580 0.2170007 1.727614e-41
## 9682 ENSG00000154545 61.40843 10.796944 3.6045524 4.601302e-41
## 1809 ENSG00000089057 1005.89933 -2.735223 0.2068275 6.323177e-41
## 10094 ENSG00000158258 1125.91097 -6.095347 0.4558747 7.334541e-41
## 16378 ENSG00000196839 229.11553 2.746156 0.2083082 1.988212e-40
## 9239 ENSG00000150760 1196.20071 -2.184512 0.1657465 5.495730e-40
## 9251 ENSG00000150907 1129.39294 -3.533868 0.2688621 7.581393e-40
## 7926 ENSG00000140577 937.43722 -2.875797 0.2199557 7.790204e-40
## 11271 ENSG00000164611 496.61379 4.635777 0.3564130 9.545710e-40
## 11873 ENSG00000167074 615.93585 -3.346519 0.2564730 1.743841e-39
## 12704 ENSG00000170653 777.55511 -1.810626 0.1394129 5.215307e-39
## 1558 ENSG00000082175 367.99463 -7.875171 0.6011838 5.722751e-39
## 8489 ENSG00000144445 421.28187 -2.167882 0.1672635 8.997096e-39
## 10100 ENSG00000158301 326.81998 -4.133615 0.3200302 9.061444e-39
## 2252 ENSG00000100320 2020.23680 -2.638264 0.2043416 9.176720e-39
## 7145 ENSG00000136111 743.08883 -4.240211 0.3274465 1.049670e-38
## 4417 ENSG00000115504 825.88386 -2.683715 0.2074528 1.110653e-38
## 18966 ENSG00000249242 299.88110 -4.701757 0.3636611 1.612567e-38
## 6320 ENSG00000131018 1294.67531 -4.585675 0.3573770 1.723938e-38
## 11795 ENSG00000166783 1382.42731 -1.865451 0.1443335 2.355278e-38
## 2831 ENSG00000103657 1028.96224 -2.041921 0.1596197 5.691320e-38
## 6335 ENSG00000131089 432.50559 -2.431705 0.1912510 1.544106e-37
## 2254 ENSG00000100324 462.88611 -1.932270 0.1523835 1.560118e-37
## 13467 ENSG00000174485 273.99563 -2.303004 0.1815890 2.447879e-37
## 15364 ENSG00000185551 3099.39881 -4.350840 0.3417257 3.205223e-37
## 7644 ENSG00000138688 1228.68677 -2.051523 0.1633682 1.434492e-36
## 399 ENSG00000020181 1145.56744 -5.417010 0.4341985 2.411756e-36
## 15709 ENSG00000187240 368.43998 -2.941907 0.2365892 4.054844e-36
## 863 ENSG00000063587 889.69004 -3.179799 0.2559716 5.350762e-36
## 5890 ENSG00000127603 4721.59003 -1.954497 0.1570793 6.035839e-36
## 9996 ENSG00000157404 334.79742 -6.776571 0.5421366 6.416311e-36
## 11981 ENSG00000167553 5217.60943 3.045495 0.2476570 6.839267e-36
## 1107 ENSG00000070366 954.70711 -2.255985 0.1818605 8.455774e-36
## 13342 ENSG00000173714 518.81452 -10.453879 0.8596890 1.146955e-35
## 16609 ENSG00000197724 860.07611 -2.069177 0.1667896 1.202119e-35
## 4147 ENSG00000113319 310.26753 -5.227584 0.4236730 1.268900e-35
## padj symbol
## 4410 8.301023e-54 IGFBP5
## 2280 2.858189e-50 ST13
## 1209 5.161453e-48 EVC
## 2216 1.465950e-47 POLDIP3
## 4194 4.209616e-44 SMAD5
## 4809 7.492388e-42 VAMP8
## 11261 1.542669e-41 GALNT10
## 1924 5.722577e-41 MAP3K20
## 14389 6.760089e-41 ZNF609
## 1090 5.449016e-40 GNB5
## 8092 6.122744e-40 <NA>
## 3841 8.109052e-40 ARPC3
## 7918 2.353624e-39 DET1
## 11027 2.300146e-38 FYCO1
## 15073 2.300146e-38 PGP
## 9682 5.743288e-38 MAGED4
## 1809 7.428246e-38 SLC23A2
## 10094 8.137673e-38 CLSTN2
## 16378 2.089820e-37 ADA
## 9239 5.487761e-37 DOCK1
## 9251 7.071735e-37 FOXO1
## 7926 7.071735e-37 CRTC3
## 11271 8.288582e-37 PTTG1
## 11873 1.451094e-36 TEF
## 12704 4.166196e-36 ATF7
## 1558 4.395733e-36 PGR
## 8489 6.319596e-36 KANSL1L
## 10100 6.319596e-36 GPRASP2
## 2252 6.319596e-36 RBFOX2
## 7145 6.987653e-36 TBC1D4
## 4417 7.155111e-36 EHBP1
## 18966 1.006393e-35 TMEM150C
## 6320 1.043296e-35 SYNE1
## 11795 1.383449e-35 MARF1
## 2831 3.247467e-35 HERC1
## 6335 8.420842e-35 ARHGEF9
## 2254 8.420842e-35 TAB1
## 13467 1.286490e-34 DENND4A
## 15364 1.641321e-34 NR2F2
## 7644 7.162062e-34 BLTP1
## 399 1.174761e-33 ADGRA2
## 15709 1.928078e-33 DYNC2H1
## 863 2.485118e-33 ZNF275
## 5890 2.739585e-33 MACF1
## 9996 2.847559e-33 KIT
## 11981 2.969283e-33 TUBA1C
## 1107 3.592984e-33 SMG6
## 13342 4.772049e-33 WFIKKN2
## 16609 4.899493e-33 PHF2
## 4147 5.068241e-33 RASGRF2final_data_Treatment <- read.csv('Final_Data_Treatment.csv')
final_data_Treatment$symbol <- mapIds(org.Hs.eg.db, keys = final_data_Ascites$X, column = 'SYMBOL', keytype = 'ENSEMBL')## 'select()' returned 1:many mapping between keys and columnsfinal_data_Treatment <- final_data_Treatment[order(final_data_Treatment$pvalue),]
head(final_data_Treatment, 50)## X baseMean log2FoldChange lfcSE stat
## 18494 ENSG00000232196 2250.560333 -5.488205 0.4186856 -13.108177
## 9507 ENSG00000152954 427.334660 -6.155589 0.4783093 -12.869474
## 18072 ENSG00000214954 473.432278 -4.931154 0.4457730 -11.062030
## 15436 ENSG00000185873 27.757096 -7.008933 0.6421647 -10.914542
## 6142 ENSG00000130032 1514.000671 -4.539955 0.4188374 -10.839420
## 9166 ENSG00000149742 25.142249 -4.488554 0.4279113 -10.489450
## 9568 ENSG00000153498 12.539244 -6.535808 0.6400188 -10.211901
## 18007 ENSG00000214306 22.749507 -4.970228 0.5133714 -9.681545
## 19229 ENSG00000256222 408.858754 -4.011969 0.4152292 -9.662057
## 19221 ENSG00000256045 4355.761910 -2.870513 0.2979582 -9.633946
## 16689 ENSG00000197978 31.963561 -21.631656 2.2510063 -9.609771
## 16200 ENSG00000196242 45.534823 -5.892379 0.6152154 -9.577750
## 8677 ENSG00000145832 137.012400 -4.511195 0.4781718 -9.434256
## 19368 ENSG00000258847 58.250413 -5.606282 0.5987555 -9.363224
## 6636 ENSG00000133063 568.796354 5.910675 0.6322847 9.348123
## 13679 ENSG00000175728 19.170085 -4.289633 0.4605539 -9.314074
## 19873 ENSG00000269831 20.630671 -4.313918 0.4640769 -9.295698
## 13666 ENSG00000175604 27.136243 -4.047076 0.4379947 -9.240012
## 14174 ENSG00000178928 11.966382 -5.648865 0.6368978 -8.869343
## 19728 ENSG00000268640 16.625329 -4.379262 0.4941821 -8.861636
## 12882 ENSG00000171471 29.545943 -2.820267 0.3230463 -8.730227
## 17033 ENSG00000203443 13.473072 -5.956947 0.6928507 -8.597735
## 15023 ENSG00000183888 40.049824 -2.817759 0.3322119 -8.481814
## 19706 ENSG00000268423 6.954908 -4.992747 0.5948706 -8.392996
## 14642 ENSG00000182177 16.331918 -2.987088 0.3574802 -8.355952
## 17336 ENSG00000204752 12.917267 -3.777498 0.4541844 -8.317102
## 16415 ENSG00000196979 24.325575 -2.246449 0.2722060 -8.252754
## 11188 ENSG00000164287 95.468586 -3.451748 0.4198622 -8.221145
## 17425 ENSG00000205147 35.643091 -2.579262 0.3141514 -8.210254
## 11540 ENSG00000165702 15.234682 -4.243058 0.5172300 -8.203427
## 11868 ENSG00000167014 14.532269 -4.798287 0.5904952 -8.125869
## 17359 ENSG00000204869 8.745445 -4.830974 0.5955747 -8.111450
## 19102 ENSG00000254440 19.383080 -3.359643 0.4147415 -8.100572
## 10528 ENSG00000161807 5.136641 -5.030228 0.6224915 -8.080798
## 8583 ENSG00000145107 75.409737 2.685922 0.3325451 8.076864
## 15462 ENSG00000185960 97.400678 -2.370067 0.2947679 -8.040452
## 8718 ENSG00000146070 217.320459 2.765592 0.3453601 8.007851
## 12427 ENSG00000169347 28.730232 -3.011684 0.3784266 -7.958436
## 17571 ENSG00000206260 28.753091 -2.604763 0.3309301 -7.871037
## 2783 ENSG00000103310 6.264883 -4.909625 0.6243926 -7.863042
## 7907 ENSG00000140505 49.809605 -2.303440 0.2968559 -7.759455
## 16417 ENSG00000196987 10.157904 -3.419154 0.4434952 -7.709561
## 18636 ENSG00000237353 18.623319 -3.096315 0.4023738 -7.695122
## 9795 ENSG00000155714 13.428164 -3.450802 0.4491391 -7.683148
## 14441 ENSG00000180777 27.742223 -2.368154 0.3089783 -7.664467
## 3262 ENSG00000106302 11.667179 -3.858102 0.5036261 -7.660647
## 19667 ENSG00000268154 17.275943 -3.590645 0.4703960 -7.633239
## 6625 ENSG00000132975 46.293578 -3.482568 0.4563831 -7.630799
## 5336 ENSG00000123342 1165.171869 1.772749 0.2325306 7.623724
## 19656 ENSG00000268046 3.907754 -4.847185 0.6384249 -7.592412
## pvalue padj symbol
## 18494 2.956193e-39 5.787339e-35 ETS2
## 9507 6.685600e-38 6.544200e-34 MTX2
## 18072 1.917078e-28 1.251021e-24 GPR149
## 15436 9.822333e-28 4.807295e-24 RAB43
## 6142 2.238858e-27 8.766023e-24 RPRM
## 9166 9.658912e-26 3.151542e-22 ZNF69
## 9568 1.753943e-24 4.905277e-21 CFAP44
## 18007 3.612144e-22 8.839368e-19 SKIC8
## 19229 4.370002e-22 9.505727e-19 CDH2
## 19221 5.747896e-22 1.125266e-18 HNRNPCL1
## 16689 7.271037e-22 1.294046e-18 LDLRAD1
## 16200 9.918170e-22 1.618067e-18 UBAC2
## 8677 3.937798e-21 5.930021e-18 ASB10
## 19368 7.734019e-21 1.081492e-17 UBXN2A
## 6636 8.921713e-21 1.164402e-17 EDN1
## 13679 1.230204e-20 1.505231e-17 CDC42EP2
## 19873 1.462436e-20 1.684124e-17 IDNK
## 13666 2.464699e-20 2.680634e-17 FBXO48
## 14174 7.357886e-19 7.581333e-16 LINC00923
## 19728 7.884895e-19 7.718129e-16 ACBD3
## 12882 2.541615e-18 2.369390e-15 <NA>
## 17033 8.130506e-18 7.235042e-15 <NA>
## 15023 2.217096e-17 1.887134e-14 <NA>
## 19706 4.739062e-17 3.865693e-14 ELOVL5
## 14642 6.490700e-17 5.082737e-14 NOVA2
## 17336 9.013995e-17 6.787192e-14 CALM2
## 16415 1.547859e-16 1.122312e-13 TRBV27
## 11188 2.015689e-16 1.409326e-13 <NA>
## 17425 2.207219e-16 1.490025e-13 MS4A12
## 11540 2.336300e-16 1.524592e-13 LTA4H
## 11868 4.441673e-16 2.804988e-13 SPARC
## 17359 5.001933e-16 3.060089e-13 CEL
## 19102 5.470135e-16 3.245116e-13 MORC4
## 10528 6.434429e-16 3.704906e-13 ZNF665
## 8583 6.645375e-16 3.717043e-13 LCN6
## 15462 8.950794e-16 4.867491e-13 KRT85
## 8718 1.167306e-15 6.176310e-13 TRIM16
## 12427 1.742268e-15 8.975892e-13 NR2C1
## 17571 3.517148e-15 1.765518e-12 ERMP1
## 2783 3.749160e-15 1.834933e-12 ZNF304
## 7907 8.529535e-15 4.072749e-12 TMEM216
## 16417 1.262510e-14 5.884801e-12 RHOQ
## 18636 1.413594e-14 6.435798e-12 ZNF526
## 9795 1.552261e-14 6.906502e-12 CIB3
## 14441 1.795755e-14 7.812331e-12 PLEKHM1
## 3262 1.849990e-14 7.873318e-12 CDS2
## 19667 2.289280e-14 9.515363e-12 SLX4IP
## 6625 2.333031e-14 9.515363e-12 HYOU1
## 5336 2.464603e-14 9.846843e-12 RGS17
## 19656 3.140049e-14 1.229455e-11 CHD5